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■ A renormalization-scale-invariant generalization of the diagonal Pade 

approximants (dPA), developed previously, is extended so that it becomes 
Q-i! renormalization-scheme-invariant as well. We do this explicitly when two 

Q_ij terms beyond the leading order (NNLO,~a^) are known in the truncated 

perturbation series (TPS). At first, the scheme dependence shows up as a de- 
pendence on the first two scheme parameters C2 and C3. Invariance under the 
change of the leading parameter C2 is achieved via a variant of the principle of 
minimal sensitivity. The subleading parameter C3 is fixed so that a scale- and 
scheme-invariant Borel transform of the resummation approximant gives the 
correct location of the leading infrared renormalon pole. The leading higher- 
twist contribution, or a part of it, is thus believed to be contained implicitly 
in the resummation. We applied the approximant to the Bjorken polarized 
sum rule (BjPSR) at Qp h = 5 and 3 GeV 2 , for the most recent data and the 

data available until 1997, respectively, and obtained af s (M%) = 0.119^0^ 
and 0.113^o'oi9, respectively. Very similar results are obtained with the Grun- 
berg's effective charge method and Stevenson's TPS principle of minimal sen- 
sitivity, if we fix C3~parameter in them by the afore-mentioned procedure. The 
central values for a^ s (M|) increase to 0.120 (0.114) when applying dPA's, 
and 0.125 (0.118) when applying NNLO TPS. 
PACS number(s): ll.10.Hi, 11.80.Fv, 12.38.Bx, 12.38.Cy 
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I. INTRODUCTION 



The problem of extracting as much information as possible, from an available QCD or 
QED truncated perturbation series (TPS) of an observable, and including this information in 
a resummed result, was the focus of several works during the last twenty years. Most of these 
resummation methods are based on the available TPS only. Some of these latter methods 
eliminate the unphysical dependence of the TPS on the renormalization scale (RScl) and 
scheme (RSch) by fixing them in the TPS itself. Among these methods are the BLM fixing 
motivated by large-nj considerations U, principle of minimal sensitivity (PMS) 0, effective 
charge method (ECH) [[||J (cf. Ref. |5| for a related method). Some of the more recent 
approaches in this direction include approaches related with the method of "commensurate 
scale relations" an approach using an analytic form of the coupling parameter M, ECH- 
related approaches [§], a method using expansions in the two-loop coupling parameter || 
expressed in terms of the Lambert function [JTTJ, methods using conformal transformations 
either for the Borel expansion parameter or for the coupling parameter ||12| . A basically 



different method consists in replacing the TPS by Pade approximants (PA's) which provide 
a resummation of the TPS such that the resummed results show weakened RScl and RSch 
dependence |13|| . In particular, the diagonal Pade approximants (dPA's) were shown to 



be particularly well motivated since they are RScl-independent in the approximation of the 
one-loop evolution of the coupling a s (Q 2 ) |14[]. An additional advantage of PA's is connected 



with the fact that they surmount the purely polynomial structure of the TPS's on which they 
are based, and thus offer a possibility of accounting for at least some of the nonperturbative 
contributions, via a strong mechanism of quasianalytic continuation implicitly contained in 
PA's. 

Recently, we proposed a generalization of the method of dPA's which achieves the exact 



perturbative RScl independence of the resummed result ||15|| . While this procedure in its 
original form was restricted to the cases where the number of available TPS terms beyond 
the leading order (LO: ~a 1 ) is odd, it was subsequently extended to the remaining cases 
where this number is even [|l(J. This would then apply to those QCD observables where 
the number of such known terms is two (NNLO,~af).0 In JT^] we also speculated on ways 



how to eliminate the leading RSch-dependence from our approximants A, and proposed 
for the NNLO case a simple way following the principle of minimal sensitivity (PMS). It 
turns out that the way proposed there does not work properly in practice since no minimum 
of the PMS equation dA/dc2 = [cf. Eq. (40) there] can be found. The dependence of 
our approximants on the RSch-parameters c 2 = /3 2 /Po and c 3 = /3 3 //3 of the original TPS 
is definitely a problem when the approximants are applied to the low-energy observables 
like the Bjorken polarized sum rule (BjPSR) at the low momentum transfer of the virtual 
photon, e.g. Q 2 ph ^3-5 GeV 2 |7 



In the present work, we address this problem. For the NNLO TPS case, we construct 
in Section II an extended version A of our approximants, in which the dependence on the 
leading RSch-parameter c 2 is successfully eliminated by application of a variant of PMS 



2 When just one such term is known (NLO), our approximants give the same result as the ECH 
method. 
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conditions dAjdSp = 0. This procedure can be extended in a straightforward way to the 
cases where more terms are known in the TPS, e.g. the NNNLO cases available now in QED, 
but we will not discuss such cases here. In Section III, we apply our approximant to the 
BjPSR at such Q^ h where three quark flavors are assumed active, e.g. Q^ h ~ 3-5 GeV 2 . While 
the approximant at this stage is an RScl-independent and C2-independent generalization of 
the diagonal Pade approximant (dPA) [2/2], it still contains C3-dependence comparable to 
that of the ECH || and TPS-PMS methods. Subsequently, we fix the value of c 3 in our, 
the ECH and the TPS-PMS approximants so that PA's of a modified (RScl- and RSch- 
independent) Borel transform of these approximants yield the correct location of the leading 
infrared (IR) renormalon pole. Thus, in the approximants we implicitly use /3-functions 
which go beyond the last perturbatively calculated order of the observable (NNLO), in 
order to incorporate the afore-mentioned nonperturbative information. In Section IV we 
then compare the values of these resummation approximants with the values for the BjPSR 
extracted from experiments, and obtain predictions for a s (M§). We also apply the TPS and 
various PA methods of resummation to these values of the BjPSR and obtain higher values 
for ot s (M|). In Section V we redo the calculations by applying PA-type of quasianalytic 
continuation for the /3-functions relevant for our, ECH, and TPS-PMS approximants. We 
further address the question of higher-twist terms. In Section VI we discuss the obtained 
numerical results for cn s (M§), and Section VII contains summary and outlook. 

A brief version containing a summarized description and application of the method can 



be found in [Tj|. In contrast to HI8] , the numerical analysis of the BjPSR in the present 
paper (Sections IV, V) uses, in addition, the most recent data of the E155 Collaboration 

m 



II. CONSTRUCTION OF C 2 INDEPENDENT APPROXIMANTS 

Let us consider a (QCD) observable S, with negligible mass effects, which is normalized 
so that its perturbative expansion takes the canonical form 

S = a (l + ria + r 2 al + r 3 al ^ ), (1) 

where ao = a4°V 7r - We suppose that this expansion is calculated within a specific RSch and 
using a specific (Euclidean) RScl Q (symbol '0' is generically attached to the RScl and 
RSch parameters in the TPS) up to NNLO, yielding as the result the TPS 

S{2] = Oo(l + r iao + r 2 a 2 ) . (2) 

Here, both a and the coefficients r\ and r 2 are RScl- and RSch-dependent. The coupling 
parameter a = a s /7t evolves under the change of the energy scale (RScl) Q, within the given 
RSch, according to the following renormalization group equation (RGE): 

<9a(ln<3 2 ; cf\ ■ ■ •) 2 {0) 2 , (o) 3 , \ / Q \ 
d ln (Q 2 ) = -Po a I 1 + c i a + 4 a +c\'a +••■), (3) 
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where /3q and c% are universal quantities (RScl- and RSch-invariant),(] whereas the remaining 
coefficients cp (j > 2) are RSch-dependent and their values can - on the other hand - be 
used to characterize the RSch. Consequently, in @ the coupling parameter a is a function 
of the RScl and RSch 



a = a(lnQo;4 0) ,4 0) , 



(4) 



The NLO and NNLO coefficients in (^) have, due to the RScl and RSch independence of S, 
the following RScl and RSch dependence: 



j~i = j~i 



(lnQl)=r 1 (lnQ 2 )+(3 ln (qI/q^ 



r\{\nQl) + anQnQl 



(o) , 



(5) 



where p 2 is RScl- and RSch-invariant. Although the physical quantity S must be indepen- 
dent of the RScl and RSch, its TPS (Q) possesses an unphysical dependence on RScl and 
RSch which manifests itself in higher order terms 



[2] 



dS 



[2] 



dS 



[2] 



d\nQ 2 



(6) 



All approximants to S which are based on TPS (§) must fulfill the Minimal Condition: 
when expanded in powers of ao to order Oq, they must reproduce TPS (0). Further, since 
the full S is RScl- and RSch-independent, the approximant should preferably share this 
property with S if it is to bring us closer to the actual value of S. The generalization of 



the diagonal Pade approximants developed in Ref. |15j possesses full RScl independence for 
massless observables. 

In its original form it is accountable only to TPS with an odd number of terms beyond the 
leading order (LO: ~a 1 ). Unfortunately, however, QCD observables have been calculated 
at most to the NNLO, i.e., at best the TPS (0) is known. Therefore, in Ref. |16[ we have 
extended the method to the cases with even numbers of terms beyond the LO, in particular 
for the TPS of the type (0). Since within the present paper we are going to apply an 
extended related procedure to these cases of Sp], we recapitulate briefly the main steps for 
treating a TPS of the generic form Sm- The trick consisted in introducing - in addition to 
S - the auxiliary observable S = S * S, which then gets the following formal canonical form: 



where 



S 

R2 



(S) 2 = a {0 + a + R 2 a 2 Q + R 3 a 3 + 

2r l5 /.', - 2 



(7) 
(8) 



rf + 2r 2 , . . . 

S is then known formally to NNNLO (~a 4 ) and the method can thus be applied, yielding 



an approximant A^{ 2 ^ to S. The corresponding approximant to S is \j A l ^ Zi which has the 



I [2/2] 



form [16 1 



4 2/2] = {« [a(lnQl;4 



.(o) jo) 

C 3 5 



a(\nQ 2 2 ;c2 



(0) to) 

5 c 3 i 



} 1/2 (= s [2] + o{4 



(9) 



Po = (11 — 2nj/3)/4, c% = (102 — 38nj/3)/(16/3o), where nj is the number of active quark flavors. 
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and it is again exactly RScl-invariant. Here, the two scales Qj (j = 1, 2) and the factor «o 
are independent of the RScl Q and determined by the identities 



MQVQl) 



1 

Wo 



Abo ' 



b\ = c\ — 2ri 



-c 2 + r (0) 

2 C l + C 2 



Ciri + 3r 2 - 2r 2 



(10) 
(11) 



If we ignore all higher than one-loop evolution effects, i.e., if we set c\ = = c% in (|10[H 
(p|) and replace the two coupling parameters in (§) by their one-loop evolved (from RScl 
Qq to Q|) counterparts, then the approximant (|S[) becomes the square root of the [2/2] 
Pade approximant of S. This follows from general considerations in [Fj^lC], but can also be 

1/2 

verified directly in this special case. The approximant [2/2] ^ preserves the RScl-invariance 
only approximately (in the one-loop RGE approximation). 

Although the RScl dependence is eliminated completely by using the approximant @, 
there remains a RSch-dependence, i.e., dependence on cf^ (j > 2). It manifests itself to 



a large degree due to db2/dc 2 ^ ^ ^2/802' = 3). In Ref. [|Tj| we speculated that the 
dependence on the leading RSch-parameter could be eliminated by imposing the PMS 



(o) 



condition of local independence (cf. Eq. (40) in |L6| ) 



dA*n f{mQf(cf )}, 



JO) 

C 2 5 



.(0) 



(0) 



,(0) 



dc\ 

where implicitly "= 0" should be understood as " 
However, expansion of this expression in powers of the coupling ao (or: any a) yields 



(12) 



,6" 



since in general this derivative is 



'0 • 



dM /2] 



(0) 



,(0) 



-lOciajj 



0(a 6 Q ) 



(13) 



This implies that the approximant (0) to S has no stationary (PMS) point with respect to 
the RSch-parameter c\ ° , since the coefficient of the leading term in the expansion of the 
derivative is constant and cannot be made equal zero by a change of the RSch. Also actual 
numerical calculations for various observables S confirm this. 

Therefore, we will modify the approximant (|9|) so that the new one will allow us to remove, 
by a PMS condition, the dependence on the leading RSch-parameter c^. This modification 
must, of course, be such that the afore-mentioned Minimal Condition is satisfied and that 
the RScl-invariance is preserved. We do this in the following way. We keep the overall 



functional structure of 



However, we replace the single set of RSch-parameters c) 



(0) 



(i) 



(j > 2), which we inherited from the TPS, by two sets of apriori arbitrary parameters Cj 

and Cj (j > 2) in the two coupling parameters, respectively, and we also admit new values 
of the reference momenta Q\ and Q\ 



U 1 



[2/2] 



OL 



a(lnQ?;c?\q?\ 



a(lnQ 



2.J2) (2) 
2) c 2 1 c 3 1 



j 1/2 



Si 



[2] 



o{4) 



(14) 
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The parameters and cf 1 will be appropriately fixed. They will turn out to be independent 

of the RSch-parameters cf^ and of the RScl Ql of the original TPS, just like the scales Q\ 

and Q\ and the parameter at will be.Q We will now require c 2 7^ c 2 > i n contrast to @ 
which led us to the problem flllf ). This requirement is not unnatural, since the forms (JD 
and ([TJ|) have Q\^Q\ and Q\^Q 2) respectively. The two new momentum scales Qj and the 



parameter a in fll4l) will be determined, in terms of c k f s (fc = 2, 3; j = 1,2), by expanding the 



two coupling parameters in power series of the original coupling ao (01) and requiring that 
the Minimal Condition be fulfilled, i.e., that the power series for A s i coincides with that of 
S (0)-© up to (and including) ~Oq. For this purpose we use the expansion for the general 
a = a(lnQ 2 ; c 2 , C3, . . .) in powers of ao = aQnQ 2 ,; c 2 , c\ , . . .) as obtained in Appendix A 
[Eqs. (|A.7|) - (|A.9|) 1, and apply it to as yet unspecified parameters Q\, Q\ and cj^ (j = l,2). 



at al : 2ri = — 



The resulting expressions, when introduced into the square of the right-hand side of (|14|) . 
yield an expansion in powers of ao- According to the Minimal Condition, it should coincide 
with (0) up to ~<2q. Comparison of the coefficients of Oq (n = 2,3,4) leads to the following 
relations: 

atQ ° : 1 = -^-^)' =► « = J^h)= p HQl/Ql) ■ (15) 

\-x 2 2 ) - c 1 (xi-x 2 ) + 5c 2 /{xi-x 2 ) , (16) 

at Oq : 2r 2 +r 2 = — — (xf — x\) + ^Ci(x 2 — x 2 ) — c 2 Q \x\ — x 2 ) 

-3(xi54 1) -x 2 5c 2 2) ) + ^5c 3 /{xi-x 2 ) , (17) 

where we have used the notations 

x^foHQj/Ql) , Scf = of - of (j = l,2), (18) 
5c 2 = c?-cf, 5c 3 = cW-4 2 ). (19) 

Eqs. (|T6| ) and ( |ITD are the two equations which determine the two scales Qi and Q 2 (<^> 
parameters X\ and x 2 ) as functions of cj^'s (k = 2, 3; j = 1, 2). In order to see that these two 
scales are independent of the original RScl (Q ) and of the original RSch (cjj. , /c > 2), we 
introduce 

% = /3 ln(Q 2 /A 2 ) = 1,2), (20) 

where A is the universal QCD scale appearing in the Stevenson equation (|A.1|) , so it is RScl- 
and RSch-invariant. After some algebra, we can rewrite Eqs. (|IE|) and (|T7|) as a system of 
equations for Xj 



4 Parameters c 2 and cf will be chosen later in the Section, by following a variant of the PMS; 
cij and c^ 1 will be set equal to each other and fixed in the next Sections. 
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2pi +Ci — (x 1 + x 2 ) + -- 



5 Co 



[X1-X2 



2p 2 + 3p 1 - 2cipi = (x 1 +xix 2 + x 2 ) - -c 1 (x 1 +x 2 ) + 3 



(~ (1) ~ (2)\ 

(X1C2 -X2C2 ; ; 
(£i-x 2 ) 



5ci 



2(x 1 -x 2 ) 



(21) 
(22) 



where p\ and p 2 ar e the usual RScl- and RSch-invariants as defined, e.g., in [cf. also 
(I)] 



Pl = f3 \n(Q 2 /A 2 )- ri 

P2 



2 , (0) 

r 2 — r l — C\T\ + c 2 



(23) 
(24) 



Therefore, Eqs. (pT|)— (p2|) show the following: If dp and c 2 z> and 5o 
and fixed, then the solutions Xj and thus the scales Qj (j = 1, 2) are independent of the RScl 
(Qo) and of the RSch (d 2 \ c^, . . .). Thus, we have 



.(2) 



■cjj 2 '' are chosen 



Q) = Q)(<£\<% ) > 8 <%) (7=1,2) 



a: 



PoHQi/Qi) 



(1) r (2) 
2 5 c 2 ) 



(25) 



Therefore, our approximant (0) will be regarded from now on as a function of only c 



parameters (k > 2; j — 1, 2): ^4^^ (c 2 , c^; c^, c^; . . .). For actually solving the equations 
for the scales Q\ and Q 2 , it is more convenient to use Eqs. (]16|)-(|17D . F° r the subsequent 
use, we rewrite them in the following form: 



V- ~ y- z o( c 2 } ) + y~j(5c 1 5c 2 - 5c 3 



16 V ; 
1 16c 2 





y-\ 



where we use the notations 



y± 



Set 



5A 

J 1 ) 



In ± In 

<5o 



.(2) 



(fc = 2,3) 



2p 2 + -cf 



3 c 



(26) 
(27) 

(28) 

(29) 
(30) 



where p 2 is given by fl24j) . Incidentally, it can be explicitly checked that in the special case 



of c« 



c (2) 

c 2 



dtp and = c 3 z> = c 3 u> Eqs. (p6l)-(^0|) and (|T^) recover the old approximant 



.(1) _ „(2) 



.(0) 



of Ref. [|ie 



The next question is how to fix parameters dp and c^ - ' (j= 1, 2). Above all, we have to 
fix the leading parameters dp's since otherwise their arbitrariness would reflect the fact that 



So) 



5 Raczka |2(J used the sum of the absolute values of terms in p 2 for a formulation of criteria 
for acceptable RScl's and RSch's in NNLO TPS. He concluded that the strong RScl and RSch 
dependence of the NNLO TPS of the BjPSR (with nj = 3) presents a serious practical problem. 
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the leading RSch-dependence (i.e., the dependence on c 2 ) has not been eliminated from 
the approximant. We do this by requiring the local independence of the approximant with 



respect to variation of c 2 and of c K 2 ' separately. This condition is a variant of the principle 
of minimal sensitivity (PMS), or a PMS-type ansatz 



OA 



[2/2] 



dA 



2/2] 



.(2) 



dA 



[2/2] 



dci s) 



dA 



[2/2] 



SC2 



d(5c 2 ) 



(31) 



Here, "=0" should be understood as "~cto" since in general these derivatives are ~<2q. These 
two equations then give us solutions for the leading parameters c 2 ^ and c 2 , once the values 



of the subleading parameters c 3 



+ 4 z; )/2 and 5c 3 = c : 



(1) J2) 

w 



have been chosen.[] 



However, using Eq. (|A.5|) and the fact that Q 2 - are independent of c 3 [cf. (|25|) 1, we can show 
the following dependence of the approximant on (at constant <5c 3 ): 



din 




(i(c 3 s) ) - (al + a\a 2 + a 1 al + al) + C(a^) ^ d^f)^ 



(32) 



(j = 1,2) and we took the index convention |ai| > \a 2 \ 



where Oj = a(ln cjf , c 3 , . 

This means that the dependence on c 3 cannot be eliminated in the considered case, not 
even by a PMS variant. In this respect, the situation is analogous to the usual TPS-PMS 
|| and the ECH |3j methods. These two methods (cf. Appendix C), while fixing RScl 
(Qo l— ► Qech = Qpms) an d c 2 RSch-parameter (c!j | — > C2 MS or cf CH ) in the original TPS 
(0), leave the value of the subleading parameter c 3 there unspecified, with the residual 
c 3 -dependence of the (TPS-) approximant 



din (Sgp) « d(c 3 )4/2 , 



(33) 



where label 'X' stands either for 'ECH' of 'TPS-PMS'. Comparing (53) and 



Jj, we see 

that the ^""-dependence of our approximant could be up to twice as strong as that of the 
TPS-PMS and ECH methods. 

Hence, varying c 3 and c 3 parameters in our approximant at this point would apparently 
not lead to any new insight. For the sake of simplicity, we choose from now on these two 
subleading parameters to be equal to each other 



c 3 



(2) _ 
C 3 = c 3 



(6c 3 = 0) 



(34) 



but we will adjust the common parameter C3 later to a physically motivated value. 

With the chosen restriction (33), the problem of finding our approximant (|i~4]) to the 
TPS (U) basically reduces to the problem of solving the system of three coupled equations 



(p6j) and (|3l|) for the three unknowns ?/_ [=/?o ^(Qi/^)] an d Sc 2 and c!j (^ c 2 and c 
For completeness, the PMS-like equations (|31"D , when £03 = = £04, are written explicitly 



Also a value of <5c4 = c 



(!) J 2 ) 



has to be chosen - see later. 



S 



in Appendix B, to the relevant order ~<2g at which we solve them - Eqs. ( p.l|) -( ]B~2"|) . 
From there and from (|26D we explicitly see that these three equations contain only the 
three unknowns (?/_, and 802) and the (known) RScl- and RSch-invariants p 2 fl24]) and 
ci =/3i//3q. Interestingly enough, these three equations do not depend on C3 (=c 3 = c!j ). 
In addition, they do not depend on any other higher order parameters c% (k > 4; j = 1, 2) 
appearing in aj = a (In Q 2 ; cjj , c 3 , c 4 , . . .), except on <5c 4 = C4 — cf^ which was taken to be 

zero in Eqs. (|B.l|) - (|B.2j ). Hence, Qj and c 2 (j = 1,2) will be functions of p 2 and ci only, 
thus explicitly RScl- and RSch-invariant. For simplicity, we want the solutions Q| and 

c 2 (j = 1) 2) to be independent of any higher order parameter (fc > 3) that possibly 
appears in our approximant, therefore we choose from now on also 5c^ (= c 4 1} -c 4 2) ) = 0. 
The solution of the mentioned three coupled equations in any specific case can be found 
numerically, e.g. by using Mathematica or some other comparable software for numerical 
iteration. Certainly we have to ensure that the program scans through a sufficiently wide 
range of the initial trial values y_ n '\ (c 2 )( in ') and (fo^)^ for iterations, in order not to miss 
any solution. The solutions which result in either \a\ ^> 1 or \6c\ -C 1 should be discarded 
since they signal numerical instabilities of the approximant [\a\ ^> 1 =>• Q\ ~ Q\ - cf. (p~5|) ] 
or are in addition physically unacceptable (|a| <C 1 =>• Q\ Q\ or Q| <C Qf). We have 
apparently two possibilities: 

• V—i c£ an d ^ c 2 are all real numbers (and thus the intial trial values as well); 

• C2 and its initial values are real; y_ and 602 and their initial values are imaginary 
numbers (c^ and are complex conjugate to each other, as are Q\ and Ql). 

In both cases, the approximant itself turns out to be real, as long as C3 is real. 

If we encounter several solutions which give different values for the approximant, we 
should choose, again within the PMS-logic, among them the solution with the smallest 
curvature with respect to and c\ ■ For such cases, we define two almost equivalent 
expressions for such curvature in Appendix B - cf. Eqs. (p.4j )- (|B.5|) . 



III. BJORKEN POLARIZED SUM RULE (BPSR): C3-FIXING 



We will now apply the described method to the case of the Bjorken polarized sum rule 
(BjPSR) pi]] . ^ is the isotriplet combination of the first moments over xej of proton and 
neutron polarized structure functions 



dx 



Bj 



9i 



ph) 



9i 



(^Bj) Qphy 



~\9a\ 

O 



S(Q 2 



ph 1 



(35) 



where p 2 = — Qph< is the momentum transfer carried by the virtual photon. The quantity 
S(Qp h ) has the canonical form ([I]). It has been calculated to the NNLO [^,^|, in the MS 
RSch and with the RScl Ql = Q^ h . The pertaining values of r x and r 2 , for those Q^ h where 

: 3), e.g. at Qp h = 3 or 5 GeV 2 , are r x = 3.5833 



three quark flavors are assumed active (n f ■ 
m and r 2 = 20.2153 [231 , so that 
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S[2](Q 2 ph ; Ql = Q 2 ph ; cf s , cf s ) = a (l + 3.5833a + 20.2153a 2 ) , (36) 
with : a = a(ln Q 2 ; of*, dp, . . .) , n f = 3 , cp = 4.471, dp = 20.99 . (37) 



The constant \qa\ appearing in (^) is known from /3-decay measurements (it is denoted 
there as \g A /gv\) 

\g A \ = 1.2670 ± 0.0035 . (38) 



Solving the coupled system of (^) and ( |B.l| )-( |B~"2l ) for the three unknowns ?/_, and <5c 2 , 



as discussed in the previous Section, results in this case in one physical solution onlyQ 
1 2 \ 

e -/3b In ^| =-1.514 (^ a = 0.3301), (39) 
c 2 s) = 3.301 , oc 2 = -3.672 = 1.465 , = 5.137 . (40) 



Parameter y + , defined in (^§), is then obtained from (pT|) . The resulting scales Qi, Q2 are 
then 0.767 GeV, 1.504 GeV (Q 2 h = 5 GeV 2 ) and 0.594 GeV, 1.165 GeV (Q 2 h = 3 GeV 2 ). 

We stress that these results are independent of the value of c 3 (0) and of c 4 and other c£' 



(k > 5; j = 1, 2) in the approximant y/As^ (IT^T) , and are independent of the choice of RScl 
Qo and RSch (cjj. , A; > 2) in the original TPS S^]- In TPS (|36|), the choice was Qo = Q P h and 



2 =c 2 vIS (=4.471). Knowing and c 2 ( j == 1 , 2 ) , for the actual evaluation of approximant 



(|I~4]) we need to assume a certain value for a fl3~7|) (at RScl Qo)- The value of a is obtained 
from (p~5| ) (a = 0.3303); the value of the coupling parameter aj = a(ln Q 2 ; c 2 , c 3 , c 4 , C5 , . . .) 
(j = l, 2) can be obtained, for example, by solving the subtracted Stevenson equation 

f3o , I1 S = i +CllI1 f T £i^ N ) + r dx <<*°+<»*> 



Qo «i V 1+C i a i/ ^° (l + c 1 x)(l+cix+4 i) x 2 +c 3 x 3 ) 

-1- Cl mf^) - r dx _ (j =i,2) . (4i) 

a U+tW y (l+ Cl i)(l+ Cl i+cfi 2 +cf¥) 

In fl4l| ) we ignored terms occ^ and higher since they are not known (c^ s is not known, 
either). Stated otherwise, we set here and in the rest of this Section: cj^ = cj^ = c^ s = for 
k > A, i.e. the /^-functions pertaining to the approximant are taken in the TPS form to the 
four-loop order. Hence, the only free parameter in the approximant VAs 2 flUD is now c 3 



[cf. condition (0)], all the other nonzero parameters (Q 2 , c^, a) have been determined and 
are c 3 - and RScl- and RSch-independent. Further, any effects due to the mass thresholds 



7 Formally, we get two solutions, but they give the same approximant, since the second solution is 
obtained from the first by Q\ <-> Q 2 and c^ 1 «-> c 2 2 \ Further, if ignoring in PMS conditions ( |B.1| )- 
(B.2) the denominators, one arrives at two additional solutions, both having c£ = (6p 2 — 7cf/4)/7; 
however, one can check that also the denominators are then zero and the derivative ( [B.lD reduces 
to 2(25c 2 — 15ciy_)ao/(3y__) which turns out to be finite and nozero. 
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(rif > 4) are ignored in (|4l|). These effects are suppressed because the difference of the two 
integrals in (|4l|) tends to cancel them. Note that the scales appearing in ( [411) (<5i~ 0.6-0.8 
GeV, Q2~l-2-1.5 GeV, Q p h = Qo~ 1-7-2.2 GeV) are all regarded to be below the threshold 
(nj = 3) I— > (nj = 4), i.e. all the active quark flavors are (almost) massless.0 

The main question appearing at this point is: Which value of C3 (= = cjj ) should 
we choose in our approximant? The two most obvious possibilities are C3 = or C3 = c^ s 
(= 20.99). The decision is far from being numerically irrelevant. If choosing for a = 
a(ln Q 2 ; cp, dp) at Q 2 ph = Q 2 = 3 GeV a typical value, e.g., a = 0.09 ap(3GeV 2 ) « 0.283, 
a^ s (M|) £30.113], we obtain the following resummed values for the BjPSR 5 

^|/ 2 ]( C3 = 0) = 0.1523 , ^ /2l ( C3 = dp) = 0.1632 . (42) 

The latter is 7.16% higher than the former. The corresponding resummed values of the ECH 
@ and TPS-PMS § are 

A^ CH (c 3 = 0) = 1535 ^ ^ECH( C3 = C MS) = o 1593 . (43) 

^p ms (c 3 = 0) = 0.1528 , ^P MS ( C3 = c !p) = 0.1588 . (44) 

The latter values (for c 3 = dp) are 3.79% (ECH) and 3.96% (TPS-PMS) higher than the 
former (for 03 = 0). Thus, the sensitivity of our approximant to the variation of C3 is in the 
considered case almost twice as large as for the ECH and TPS-PMS methods, as anticipated 
in fl32"l)-(|33|) in the previous Section. The true value of C3 in ^4| CH should be equal to p 3 , 
i.e. the third RScl- and RSch-invariant of the BjPSR, but this value is not exactly known 
because the N 3 LO coefficient r 3 in the perturbative expansion of the BjPSR is not known yet. 
The stronger C3-sensitivity should not be regarded as a negative feature of our approximant, 
but rather within the following context: 

Our approximant contains two (RScl-invariant) energy scales Qi, Q 2 - Since the consid- 
ered observable is close to the nonperturbative sector (Q p h < 2.5 GeV), the relevant scales 
Qj (~<5 P h) are low: Qi~ 0.6-0.8 GeV and Q 2 ~ 1-2-1.5 GeV. Thus the relevant coupling pa- 
rameters a,j = a (In Q 2 ; c£\ c 3 ) are large: a\ ~0.19 and a 2 ~0.11 (when c 3 is set equal c^ s and 
a = 0.09, Ql = Q 2 h = 3 GeV 2 ). Therefore, the contribution of the c 3 -term on the right-hand 
side of the integrated RGE (|4T| ) [•<=>■ differential RGE (0)] at such energy scales is not negligi- 
ble. This feature, to a somewhat lesser degree, can also be seen in the ECH and TPS-PMS 
approaches, where Qech (=Qpms) ~0.8 GeV and aECH = a(lnQ| CH ; cf CH , c 3 ) ~ 0.16 (when 
C3 is set equal to c^ s , and a = 0.09, Ql = Qp h = 3 GeV 2 ). The significant c 3 -dependence 
of all these approximants, at fixed do, reflects the fact that the coupling parameters a(Qj) 
appearing in the approximants are not small and that consequently the considered observ- 
able is in the low-energy regime. The values of Pade approximants (PA's), when applied 
to NNLO TPS of an observable (e.g., BjPSR), are also C3-dependent. However, the lat- 
ter C3-dependence, in contrast to that in the afore-mentioned approximants, is not playing 



8 In the whole paper, we ignore any quark mass effects, except later in the evolution o^ IS (Q 2 ,i 1 ) l— * 
a^ s (M^) where the quark mass thresholds are significant and accounted for. 
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a highlighted role, since the PA's depend in addition on the leading RSch-parameter C2 
(■<=> c 2 ) and even on the RScl Qq. 

The above considerations, however, do not address the important problem presented by 
([42|): Which value of parameter C3 should we use in our approximant? 

We note that C3 characterizes the N 3 LO term in the corresponding /^-function (^), and 
the information on its value in a considered approximant cannot be obtained from the NNLO 
TPS on which the approximant is based. To determine the optimal value of C3 in an approxi- 
mant (our, ECH, or TPS-PMS), an important known piece of (nonperturbative) information 
beyond the NNLO TPS should be incorporated into the approximant. There are at least 
two natural candidates for this: the location of the leading infrared (IRi) and ultraviolet 
(UVi) renormalon poles, i.e., the positive and negative poles of the Borel transform Bg(z) of 
the observable closest to the origin (- for a review on renormalons, see ||25|| ). In the case of 



the BjPSR, these two locations are known from large-/3o (large-n/) considerations [|26"| , |27 

Zpale = l/MlRl),Zpole = -l//3o(UVl)- 

Which of the two leading renormalons is numerically more important in the BjPSR 
case? In the simple Borel transform of the BjPSR, with MS RSch and RScl Qq = Q p h 
(n/ = 3), the ratio of the residues of the IRi and UVi poles in the large-/?o approximation is 
2 exp(10/3) ~ 56 ^> 1 • This would suggest strong numerical dominance of the IRi over 

UVi. However, when using there the V-scheme i.e. MS with RScl Qq = Q p h exp(— 5/6) 
(~Qech), this ratio goes down to 2. This would suggest that the UVi (vis-a-vis IRi) is 
not entirely negligible. The authors of Ref. p7 used the 't Hooft RSch and varied the RScl 
in such an approach (large-/?o, simple Borel transform, principal value prescription), and 
their Fig. 2 for the BjPSR at Q^ h = 2.5 GeV 2 suggests that IR renormalon contributions 
to S(Qp h ) are 3-4 times larger than those of the UV renormalons. The relative strength 
of the UV vs IR renormalon contributions, in the RScl- and RSch-noninvariant approach 
with simple Borel transform, appears to depend in practice on the choice of the RScl and 
RSch. Incidentally, a consideration of the status of the renormalon contributions and of their 
scheme-dependence was made in Ref. |p8[| . The question of the relative suppression of the 
(leading) UV renormalon contributions in RScl- and RSch-invariant resummations would 
deserve a further study. An additional uncertainty resides in the fact that the residues, in 
contrast to the renormalon pole locations, change and thus attain unknown values when we 
go beyond the large-/?o approximation. For the UV renormalons, this uncertainty shows up 
in an especially acute form . 

The afore-mentioned works, however, suggest strongly that, in the BjPSR case S(Q"i h = 
3-5GeV 2 ), we should preferably fix the value of C3 in our, ECH, and TPS-PMS resummation 
approximants by using IRi (z pole = l//3 ) and not UVi (z po \ c = —l/fl ) information. The 
IRi pole location can be transcribed as y po \ e = 2, where y = 2j3$z. This corresponds to 
possible renormalon-ambiguity contributions ~1/Qph to the BjPSR observable which are 
nonp ert urbative . 

We will present now an algorithm for adjusting approximately the value of C3 in our 
approximant for the NNLO TPS (fj). Briefly, it consists of the requirement that C3 must be 
adjusted in such a way that the Borel transform of the approximant has the correct known 
location of the lowest positive pole, where the latter location is obtained by construction of 
Pade approximants (PA's) of the Borel transform. 

A first idea would be to use simple Borel transforms. We would first expand our ap- 
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proximant (with a general yet unspecified c 3 ) in power series of a coupling parameter, say 

ao = a(lnQo; c 2 ) c 3 )•••)' U P ^° a certain order ~ a 3 +1 (j > 3), then obtain from this pre- 
dicted S\j] TPS the corresponding ByAz) TPS (up to of the simple Borel transform as 
schematically described by 



A l ~ /2] (a ; c 3 ) = 5g-(a ;c 3 ) = a 1 + na + r 2 a^ + rf'(<%)al H hrj r '(c 3 )a£ , (45) 

=> B -(,;c 3 ) = l + ^ + ^ + ll^^ + ... + !fM-. (46) 



The (approximate) pole structure of the simple Borel transform can be investigated by 
constructing various PA's of its TPS fl46|). The requirement that the lowest positive pole 
be at y(= 2(3$z) = 2.0 would then give us predictions for c 3 . However, this approach is in 
practice seriously hampered, because coefficients Tk/k\ of the simple Borel transform £>(z; c 3 ) 
depend very much on the choice of the RScl (Qq) and RSch (<4 , <4 , . . .). For example, if 
expanding our approximant JA^a^c^) up to ~ in an RSch with = d^ s and an 



arbitrary c 3 , and keeping the RScl unchanged (=<3ph)> we reproduce in the BjPSR case 
the first two coefficients t\ and r 2 of (|36|), while the predicted r 3 in this RSch is 



f- = 125.790... - ^- + c 3 . (47) 



-.3 
2 

The PA's [2/1] or [1/2] of the corresponding simple Borel transform TPS B?^(z) would 
therefore be functions of (— c 3 /2 + c 3 ), and the requirement 2/ po ic = 2.0 would at this level 
give us only a prediction for + c 3 ), not for c 3 itself.0 For example, working with 

B^(z) in the RSch with c 3 = results in a prediction for c 3 that is by about 10.5 lower 

than the one when 4° } = dp(^ 21) is used. If using the ECH a ECH (c 3 ) @ or TPS-PMS 
Spms(c 3 ) [0 approximants instead of our approximant (where c 3 is the arbitrary subleading 
parameter used in oech and apMS - cf. Appendix C), the corresponding prediction with 
Qo = Qph is: r 3 pr ' ) = 129.8998. ..+(— c 3 °^+c 3 )/2. Hence, also in the case of these approximants 
we end up with the same kind of problem of strong RSch-dependence (c^-dependence) of 
the predicted values of c 3 . 

Therefore, we will use a variant of the RScl- and RSch-independent Borel transform 



B(z) introduced by Grunberg |30|, who in turn introduced it on the basis of the modified 



Borel transform of the authors of Ref. Ml 



rod f , 

S(Ql h ) = j Q dzexp [- Pl (Q 2 ph )z\ B s {z) . (48) 



Here, p\ is the first Stevenson's RScl and RSch invariant (|23|) of the observable S: 



Pi(Q 2 ph ) = -n(Q 2 ph /Ql)+(3o\n^=p \n^ , (49) 



9 The [1/1] PA of the simple Borel transform is independent of c 3 and of C3 . In the BjPSR case, 
in the MS RSch and at RScl Qq = 3 or 5GeV 2 , where n/ = 3, it predicts y po i c «1.6. 
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where A is the universal scale appearing in the Stevenson equation (|A.1|) , while A is a scale 
which depends on the choice of the observable S. But A is independent of RScl Qo and of 
RSch and even of the process momentum Q p ^. We note that pi(Qp h ) is, up to a constant c 
(the latter is irrelevant for the position of the poles of B$), equal to 1 / a^ loop \Qp h ) . Thus, 
B,s{z) of (|48|) reduces to the simple Borel transform, up to a factor exp(cz), if higher than 
one-loop effects are ignored. The positions of the poles of B s (z) of fl4"gp are the same as 
those of the simple Borel transform. The coefficients of the power expansion of Bs{z) of 
([|8|) are RScl- and RSch-invariant, in contrast to the case of the simple Borel transform. 
These invariant coefficients can be related with coefficients r n of S with relative ease in a 
specific RSch Ck = c\ (fc = 2,3,4, . . .), while keeping the RScl Ql unchanged 



N (fn-Cif n _i, 



B s (z) = (c lZ r z exp(-nz) E ^n+l+Xz f = ^' B ^ z ) " ( 50 ) 



Here, r n is the coefficient at a n+1 in the expansion of S in powers of d = a{\nQl; cf, cf, c\ 
and by definition f_i = 0, fo = 1. In (|50|) , we introduced the modified RScl- and RSch- 
invariant Borel transform B s (z), by extracting the factor [c\z) CxZ whose behavior at z^O 
may be problematic for PA's to deal with.0 The obtained coefficients of the power expansion 
of B s (z) are explicitly RScl- and RSch-invariant, depending only on the invariants pj (j > 
2), on Ci and on some universal constants. 

We will now calculate the invariant Borel transform B of our approximant. The 
coefficients as predicted by our approximant (14) \fAs^{c^) are functions of the only 
unknown C3 = 7^(03), k > 3], They can be obtained as coefficients of the power expansion 



of y/As2(cz) in powers of a. Looking back at the form (|i~4]) of our approximant, such a 
power expansion requires first the separate expansions of a\ = a(Ql;c2,c 3 ,0,...) and of 



a 2 = a(Q|; c| , C3, 0, . . .) in powers of a. The latter expansions can be read off Eq. ( |A.7| ), up 
to ~a 5 (there: a t— > a\ or 02; and a 1— > a.) In fact, we carried out the latter expansion up to 
~a 8 (with the help of Mathematica), which allowed us to write the approximant -\A4s 2 ( c 3) 
up to ~ a 7 . This in turn leads us to obtain the invariant Borel transform B^j(z) up to 
~2 6 , according to (|50|), and allows us to construct PA's of the Borel transform of as high 
order as [3/3], [2/4], [5/1]. The coefficients starting at z 3 are predictions of the approximant 
and are c 3 -dependent: B s (z) = l + b 1 z + b 2 z 2 + b 3 (c 3 )z 3 + ■ with 6x^—0.7516, 5 2 ~ 0.4209, 
^3( c 3) ~ (— 2.664 + 0.1667C3), etc. Construction of various PA's of that Borel transform and 
requirement that the smallest positive pole equal 2/ po ie(=2/3o%>oie)= 2.0 gives us predictions 
for C3 which are listed for the described case in the second column of Table Q. In the column 
we included values of C3 with small nonzero imaginary parts and Re(c3) ~ 10-12, since for 



10 Grunberg's [|30| Borel transform £>( Gr ') was chosen by convention as: B^ Gr, \z) = T(l + 
c\z) exp(ciz)B(z). In this way, B( Gr ^(z) ~ B(z) yj2-KC\z when z — > 00, and the coefficients of 
the power expansion of i?( Gr ) in z depend only on the RScl- and RSch-invariants pj (no de- 
pendence on ci and on T-function-related constants). We decided not to follow this convention, 
primarily since T(l + ciz) introduces spurious poles on the negative axis, the one closest to the 
origin being y(= 2[3qz)k> — 2.53. Such spurious poles not far away from the origin can significantly 
limit the PA's ability to locate correctly the leading IR renormalon pole (y po ie RJ 2.). 
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such values the PA-g's and the TPS of B are almost real, with imaginary parts less than one 
per cent of the real part for y < 1.9. In the latter cases the real part of c 3 may be regarded 
as the suggested value. The actual value of c 3 must be exactly real, but since a specific PA 
predicts only an approximate value of c 3 , this latter value is not necessarily exactly real. 
We did not include some other solutions which differ a lot from those given in the column. 



PA^ 


c 3 (V-Asp) 


c 3 (ECH) 


c 3 (TPS-PMS) 


[2/1] 


21.7 


35.1 


35.1 


[3/1] 


13.7 


19.5 


19.0 


[4/1] 


11.1 


14.4 


13.1 


[5/1] 


9.3 


11.2 


8.8 


[1/2] 


12.8 


17.3 


17.3 


[2/2] 


12.4 


16.9 


16.2 


[3/2] 


11.7±3.4i 


15.8±6.4i 


15.4±7.4i 


[4/2] 


10.3±2.8i 


12.9±5.1i 


11.6±6.8i 


[1/3] 


12.4 


16.9 


16.2 


[2/3] 


12.9 


17.4 


18.3±0.8i 


[3/3] 


10.6±2.9i 


13.6±5.5i 


12.6±7.0i 


average 


ps 12.5 


Ps 17. 


~ 16. 



TABLE I. Predictions for C3 in our, ECH and TPS-PMS approximants, using various PA's of the 
invariant Borel transform B{z) of the approximants and demanding that the lowest positive pole be at 
z P oic = l/A) (= 4/9). The higher order parameters (k > 4, j = l,2) in our approximant, and (k > 4) 
in ECH and TPS-PMS, were all set equal to zero. 



Predictions of PAg's of the intermediate orders ([3/1], [4/1], [2/2], [3/2], [1/3], [2/3]) give 
us the average value c 3 ~ 12.5 which we will adopt. The prediction by PA [2/1] differs 
from most of the other predictions, apparently because [2/1] is of low order. Predictions by 
the highest PA's ([5/1], [4/2], [3/3]) also differ from the average. The reason for this lies 
probably in the fact that these PA's contain information on many higher order coefficients 
f n (^ = 3,4,5,6) which are not contained in the TPS Sm on which the approximant y/As? 
is based. In addition, these high order PA's are implicitly dependent on the high order 
parameters and (k = 4, 5, 6, 7) which were here simply set equal to zero (we will come 
back to this point later in Section V). 

Completely analogous considerations produce the values of c 3 parameter in the ECH and 
TPS-PMS approximants. For details on the ECH and TPS-PMS methods, when applied 
to the NNLO TPS Sm (0), we refer to Appendix C. Also in this case, we make for the 
corresponding /3-functions the simple TPS choice: ECH RSch= (p 2 , c 3 , 0, ...); TPS-PMS 
RSch= (3p 2 /2, c 3 , 0, ...). The obtained predictions for c 3 for these approximants are included 
in Table Q. Again, PA [2/1] and the highest order PA's appear to give unreliable predictions. 
On the basis of the predictions of PAg's of intermediate order, we will adopt the value 
c 3 = 17 for the ECH case, and c 3 = 16 for the TPS-PMS case. The actual values of c 3 must 
be exactly real. 
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In fact, we can apply this method of determining the C3-parameter of our approximant 
(and of ECH and TPS-PMS approximants) to any QCD observable given at the NNLO 
and whose leading IR renormalon pole is known via large-/? considerations. The method, 
however, is well motivated only if there are indications that the leading IR renormalon 
contributions to the observable are larger than those of the leading UV renormalon. We 
wish to stress that our approximant, as well as the ECH and TPS-PMS approximants, 
are completely independent of the original choice of the RScl and RSch in the TPS of the 
observable, because the parameter C3 is RScl- and RSch-invariant since it is determined by 
using the RScl- and RSch-invariant Borel transform B(z). 

A few remarks about the multiplicity of the discussed IRi pole are in order. The simple 
Borel transform J2 r kZ k /k\ of S(Qp h ) behaves near z po \ e (=1/ (3q) as ~l/(z po ic— z) K where the 
multiplicity is p2| , p4| , p9| k= 1 + (/3i//3o)%, ic + (l/Po), and 7 is the one-loop anomalous di- 
mension of the corresponding two-dimensional operator appearing in the Operator Product 
Expansion for 5* (usually 7 > 0). On the other hand, the RScl- and RSch-invariant Borel 
transform (|50|) behaves near z po \ e with the simpler pole multiplicity |3l| k= / Po)- To 



our knowledge, the anomalous dimension 7 is not known in this case. However, in the case 
of the Adler function (logarithmic derivative of the correlation function of quark current op- 
erators), the one-loop anomalous dimension of the four-dimensional operator corresponding 
to the lowest IR renormalon pole there (z po \ e = 2//?o) is known [32,32] to be 7 = 0. If 7 = 
also in the BjPSR case, then the RScl- and RSch-invariant Borel transform (^8|) -(p0|) has 
k = 1, i.e., the leading IR renormalon pole is a simple pole, in contrast to the simple Borel 
transform where k is noninteger. In such a case, we may have an additional incentive to use, 
instead of the simple Borel transform, the invariant Borel transform (f48|)-(|50|) in conjunc- 
tion with the afore-described PA's of Table |. Namely, PA's are very good at discerning the 
location of a pole if such a pole is simple, and are somewhat less successful in this job if the 
pole is multiple or with noninteger multiplicity. 



IV. BJPSR: PREDICTIONS FOR THE COUPLING PARAMETER 



Now that we have fixed the values of the C3-parameter in the approximants VAs 2 ^; C3), 
ECH and TPS-PMS, the only adjustable parameter in them is the numerical value of cto = 
'(Qp h )/7r, at such Qp h where three flavors are assumed active, e.g. at Q ph = 3 or 5 GeV 2 



This ao can be obtained by requiring that it should reproduce the experimental values for 

S^Jof 



The questions connected with the extraction of the values of the BjPSR 
integral fl35|) from the measured polarized structure functions are at present not quite settled. 
One source of the uncertainty arises from the fact that these structure functions have not 
been measured at small values of xej and that, therefore, a theoretical extrapolation to 



such small xej-values is needed. The authors of ||35| , |36[| used the small-xej extrapolation as 
suggested by the Regge theory, the assumption made also by various experimentalist groups 
before 1997. The values thus obtained by |35| , |36| , on the basis of measurements at SLAC 
and CERN before 1997, are 



(Regge) 



1. 

r\9a\ 
6 



1 - 5(Q^ h = 3GeV 2 



0.164 ±0.011 



(51) 



On the other hand, the authors of |37| used a small-a^j extrapolation based on the NLO 
version of the DGLAP equations (pQCD) as opposed to the Regge extrapolation (cf. also 
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Ref. |[38|| ). This leads to higher values and larger uncertainties of the BjPSR integral. The 
values extracted in this way by |37j (their Table 4), based on SLAC data, are 



(!!)■■ q\9a\ 



S(Q 2 



P h- 



:3GeV 2 



0.177 ±0.018 



(52) 



Furthermore, most of the experimentalist groups have adopted, since 1997, similar NLO 
pQCD approaches to the small-XBj extrapolation, e.g. SMC Collaboration |39|] at CERN, 
E154 H and E155 fl|] Collaborat ions at SLAC. The most recent and updated measure- 
ments of the polarized structure functions are those of Ref. [[HJ. Their combined value of 
the BjPSR-integral at <2 2 h = 5 GeV 2 is 



(-0 



r\9a\ 

6 



1-S{C$L = 5GeV 



0.176 ±0.008 



(53) 



Apart from the problem of the small-XBj extrapolation, there is a problem of accounting for 
nuclear effects. Since the extraction of the gi structure function is based on the measure- 
ments of the structure functions of the deuteron and 3 He, nuclear effects have to be taken 
into consideration. The (multiplicative) effects due to the nuclear wavefunction have been 
taken into account in (|5^) and (|52|). However, recently the authors of [[IT] argued that addi- 
tional nuclear effects, originating from spin-one isosinglet 6-quark clusters in deuteron and 
helium (which include the shadowing, EMC and Fermi motion effects), affect the extracted 
values of the neutron structure function g± in such a way that the value of the BjPSR 
integral increases by about 10%. This would then change the E155 values of to 



(I') 



6 



1-S{CJL =5GeV 2 ) 



0.193 ±0.009 



(54) 



The values of the (H), at Q 2 h = 3GeV 2 would be increased to about 0.195 ± 0.020. We will 
not consider this case II' and case I' ( |54"D for the time being, but will briefly return to them 
in Section VI. 

In the following we will extract the values of a^ 1 (Qph) from the BjPSR-integral values 
(|53|) and (|52"D, and will simply denote the corresponding cases as I and II, respectively. 

If we insert the value (BBf) for \g^\ into fl5B|) and (|5"2"D, we obtain 



(!)■■ 
(II): 



S(Q: 
S(Q 2 



2 

ph - 



ph 



:5 GeV 2 ) 
:3 GeV 2 ) 



0.167 ±0.038 
0.162 ±0.085 



(55) 
(56) 



The present small uncertainty in the value of 
uncertainties of S(Q 2 h ) in (p5|) -(^6l). 



practically does not contribute to the 



a 



MS 



Our approximant gives, for example, for = a(ln3GeV ; C2 , C3 , 
((/- 



) = 0.09 

3GeV 2 ) ~ 0.283] the value 0.1585, which is not far from the middle values in 
(^5|)-(|56|). Varying a in our approximant (with c 3 = 12.5) in such a way that the middle 
and the end-point values of the right-hand side of ( [55"|) or ([56]) are reproduced then results 
in the following predictions for a s (in MS RSch): 



a 



MS 



(Q 2 = 5 GeV 2 



0.2894+ - 0238 



'-0.0345 



(/) 



a 



MS 



(Q 2 = 3 GeV 2 



0.2855±g;°«2 W 



(57) 
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We then evolved these predicted values via four-loop RGE (|) to Q 2 = M|, using the values 



of the four-loop coefficient c^{nf) in the MS RSch [42] and the corresponding three-loop 



matching conditions ||3 for the flavor thresholds. We used the matching at /x(%) = Km q (nf) 



with the choice k = 2 } where m q (rif) is the running quark mass m q [m q ) of the n/th flavor and 
n{rif) is defined as the scale above which rif flavors are active.^ The resulting predictions 
for a s (M|) are 

af (Ml) = 0.1196±« (/) ; 0.1135+™ (II) . (58) 

In Table [|, we give the values of as predicted from the BjPSR data (|55|) and (|56l ) by 
our approximant (with 03 = 12.5), by the ECH (with 03 = 17), and by the TPS-PMS (with 
C3 = 16). For comparison, we include predictions of these three approximants when C3 in 
them is set equal to zero, i.e., for the case when the location of the leading IR renormalon 
(IRi) pole in these approximants is not correct. Given are always three predictions for a s , 



approximant 


a,(5GeV 2 ): (I) a s (3GeV 2 ): (II) 


a s (M 2 ): (I) a s {Ml): (II) 


NNLO TPS 


n o 987 +0.04ti5 0901+0.0989 
U - JZ °' -0.0530 u - ozzi -0.1341 


n ior;o+0.0055 n 11 SO+0.0095 
U.lZOZ_ 007g U.llOO_ 232 


N 3 LO TPS (r 3 = 128.05) 


n 01 91 +0.0393 n o nfi r+0.0823 
u.jizi_ 0464 u.juoo_ 1215 


n19 on+0.0050 niiRq+ UUU8y 

U.1ZOU_ 00073 U.llDO_o.0219 


[l/2] s (NNLO) 


0.3054+-^ 0.3003+-^ 


nl99n +0.0o46 n -, -, C-K+0.U079 


[2/l] s (NNLO) 


0.3006+-^ 0.2959+-- 


nl91 o+0.0044 nluq +0.U07S 


^/[2/2} S 2 (NNLO) 


0.2937+-- 0.2895^ 


ni2n o+o.uo39 n lun +u.uuee 


[2/2) s (N 3 LO, r 3 = 128.05) 


0.2944™ 0.2901+™ 


nl o n4 +0.U040 q -| 141 +0.U069 


TPS-PMS (NNLO, c 3 = 16.) 


0-29071-^ 0.2867+^ 1035 


0-H9818S 0.1136l-. 0197 


ECH (NNLO, c 3 = 17.) 


u.2»y»_ 34 8 u.2»by_ 01028 


n 1 1 qfi+ u - UU37 n 1 1 qc;+u.oo6o 
u.iiyo_ 00059 u.iioo_ 00196 


^A [ y ] (NNLO, c 3 = 12.5) 


0.28941- 2 ! 0.2855±g:2g2 


1 1 O«+0.0035 q 1 1 oc+0.0058 


TPS-PMS (NNLO, c 3 = 0.) 


n 9Qt;7+0.U296 r\ oqi o+'- 
u - zy,J ' -0.0380 u - zyio -0.1077 


1206 +lU)U42 1143+- 

U.1ZUD_ 0064 U.114O_ 0203 


ECH (NNLO, c 3 = 0.) 


0.2947™ 0.2904™ 


i904+ u - UU3y 1142+u-uotifci 

U.lZU4:_o 0062 u - ii4z -0.0202 


(NNLO, c 3 = 0.) 


0.2960™ 0.291618:5^ 


120fi+ - 0040 ni143+ 00067 

U.lZUO_o.o063 u - ii4,; >-0.0202 



TABLE II. Predictions for a^ s , derived from various resummation approximants to the BjPSR at 
Q 2 hoton = 5GeV 2 , 3 GeV 2 . Predictions for the case I (|5^) and II ( |52| ) are given in parallel. 



corresponding to the three values of 5* (p5|) for case I, and (p6|) for case II. In addition, 
predictions of the following approximants are included in Table |lj TPS Sp] ([36]) (NNLO 
TPS); TPS S [3] with r 3 = 128.05 (N 3 LO TPS); off-diagonal Pade approximants (PA's) [1/2] <? 
and [2/1)5, both based solely on the NNLO TPS Sp] (|36|); square root of the diagonal PA 



11 If increasing k from 1.8 to 3 in case I, the predictions for the central, upper, lower values of 
a s (M§ ) decrease by 0.12%, 0.15%, 0.09%, respectively; increasing k from 1.5 to 3 in case II, the 
respective numbers are 0.12%, 0.17%, 0.03%. We assumed m c (m c ) = 1.25 GeV and 7Bb(mf,) = 4.25 
GeV. 
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(dPA) [2/2] S 2, which is based solely on the NNLO TPS ©; [2/2] s is the dPA constructed 
on the basis of the N 3 LO TPS S [3] with r 3 = 128.05. For [2/2} s and N 3 LO TPS we chose 
the latter value of r 3 (in MS, at RScl Ql = <3 2 h , rif = 3) because then the [1/2] PA of the 
invariant Borel transform B s fl50|) predicts the IRi pole y po i e = 2.0. We wrote in Table [II 
numbers with four digits in order to facilitate a clearer comparison of predictions of various 

methods. 

From Table [TT] we see that the values of a^ s (M|) predicted by various approximants 
differ significantly from each other. Addition of the N 3 LO term in the TPS decreases the 
central value of ap(Mf ) by 0.0022 (0.0020 in case II), and application of the NNLO dPA 

1 /2 

approximant [2/2] s2 decreases this value by a further 0.0027 (0.0023). Our approximant 
V As 2 (c3 = 12. 5), which is an RScl- and RSch-invariant extension of the method of the dPA 
[2/2]J/ 2 2 , decreases the central o-f^Mf ) by a further amount of 0.0007 (0.0005). Predictions 
of the ECH and TPS-PMS methods are very close to those of our method if the value of 
C3 in them is adjusted in the afore-described way. However, predictions of these two and of 
our method increase and come closer to the predictions of the NNLO dPA once we simply 
set in these approximants C3 = 0, thus abandoning the requirement of the correct location 
of the IRi pole. The predictions of the N 3 LO dPA [2/2) s are almost identical with those of 
the NNLO dPA. All the PA resummations were carried out with the RScl Qo = Qp\ 1 (n/ = 3) 
and in MS RSch, and their predictions would change somewhat if the RScl and RSch were 
changed - in contrast to the presented predictions of y/As?, ECH and TPS-PMS. 

We wish to point out that the a^ s -predictions for the case II ( |52"D were already presented 
in the short version ||18|| . However, they were somewhat lower there [the central values of 



a 



MS 



(M§) were lower by about 0.0009-0.0011] - because the value of the /?-decay parameter 
\g&\ there was taken from the Particle Data Book of 1994 \g^\ = 1.257(±0.2%) (used also in 
PTfl), while the value used here fl38|) is the updated value based on [pi] . 

In Fig. |l](a) we present various approximants for S(Q^ h ) as functions of a^ s (Qph) ( n f — 3, 
e.g. Qp h = 3 or 5GeV 2 ), and in Fig. p](b) the approximants for S^SGeV 2 ) as functions of 

a^ s (M|). There is one peculiarity of the (NNLO) TPS-PMS method, as seen also in Figs. [I] 
- for high values of observable S this method does not give solutions. This is so because the 
polynomial form of the (NNLO) TPS-PMS S FMS [see Eq. ( |C.4|) 1 is bounded from above by 
Sm&x = (2/3) 3 / 2 p2" 1 ^ 2 which, in the considered case (p2 = 5.476), is equal to 0.233 which is 
below Smax = 0.247 in case II (cf. Appendix C for more details). This is also indicated in 
Table 0. 

We wish to emphasize one aspect that makes the approximant VAs 2 conceptually quite 
different from the dPA [2/2] 5. Although both approximants incorporate information about 
the location of the IRi pole (y po \ e = 2), they do it in two very different ways. The dPA 
[2/2]5 is constructed on the basis of the N 3 LO TPS with r% = 128.05, where only this latter 
coefficient contains approximate information on the pole's location. So this dPA is a pure 
N 3 LO-construction and is RSch- and even RScl-dependent (weakly). The approximant 
y/As2 is constructed on the basis of the NNLO TPS. It is a RScl- and c^-independent 
NNLO-construction, and the correct IRi pole location is obtained by the adjustment of the 
C3-parameter within the approximant. As argued previously [cf. 2nd paragraph after ([44])], 
the c 3 -dependence in \fA~s^{c 3 ) is closely related with the sensitivity of the approximant 
to the details of the RGE evolution, and the latter details are the more important the 
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more nonperturbative the observable is. So it seems very natural that it is the intrinsic 
c 3 -parameter in y/As^{c^) that parametrizes the (nonperturbative) IRi pole location, and 
at the same time it makes the approximant fully RSch-independent. The same is true for 
the ECH and the TPS-PMS approximants. 

On the other hand, it would be an ambiguous approach to implement this kind of C3- 

1 /2 

fixing in the NNLO PA methods ([l/2]s, [2/1] g, [2/2] s ' 2 ) - because these resummations 
depend in addition on the leading RSch-parameter c 2 c^) and even on the RScl Ql- 
Therefore, it may not be so surprising that the results of our method, ECH, and TPS-PMS, 
with the mentioned C3-fixing, all give predictions that are clustered closely together and are 
significantly distanced from the predictions of (d)PA's. 

There is another theoretical aspect which indicates that the predictions of the (NNLO) 
approximant VAs 2 should in general be better than those of the (NNLO) dPA [2/2]g2 2 - 
Namely, the latter dPA is just a one-loop approximation to our approximant. More specif- 
ically, dPA [2/2]^ 2 is like ansatz (|H|), but each ctj = a(lnQ|; c\ P, c 3 , . . .) is replaced by 
the coupling parameter a( 1_1 -)(lnQ 2 ) evolved from the RScl Ql to a by the one-loop 
RGE in the original (MS) RSch. This follows from considerations in [[Tj^,[nj, and can also 
be checked directly as indicated in the parag raph after Eqs. The dPA [2/2] ^ 2 

possesses residual RScl-dependence, and RSch-dependence, the unphysical properties not 
shared by the true (unknown) sum. The approximant VAs 2 ; however, possesses RScl- and 
RSch-independence, and is thus better suited to bring us closer to the true sum. 

On the other hand, when compared with the structure of the ECH and TPS-PMS 
approximants, a/ As 2 possesses a theoretically favorable "PA-type" feature that the other 
two methods don't have: It represents an efficient quasianalytic continuation of the NNLO 
TPS S[2] from the perturbative (small-a) to the nonperturbative (large-a) regime. This is so 
because \/«4s 2 is related with the mentioned dPA method [2/2] ^ (see above). The ECH and 
the TPS-PMS approximants don't possess this strong type of mechanism of quasianalytic 
continuation, because they don't go beyond the polynomial TPS structure of the original 
TPS S[2\. These two approximants do possess, however, a weaker type of quasianalytic 
continuation mechanism, provided by the RGE-evolution of the coupling parameter a itself. 
In the one-loop limit, this would amount to the [1/1] PA-type quasianalytic continuation 
mechanism for a itself, which may explain why especially the ECH method appears to do 
well even in the deep nonperturbative regime (where S has large values). 

The possibility to adjust the value of the N 3 LO coefficient r 3 of fl3"6|) by the IRi pole 
requirement y po \ e (= 2(3 z) = 2 in the BjPSR was suggested by the authors of Ref. fl36| . 
They chose r 3 (at RScl Ql = Q^ h and in MS RSch) approximately so that the PA [2/1] 
of the simple Borel transform of that TPS gave y po \e ~ 2. In fact, they chose r 3 = 130., 
which would correspond to their 2/ po ie ~ 2.10, and then resummed the obtained N 3 LO TPS 
for S(Q 2 ph = 3GeV 2 ) by the [2/2] dPA. However, as we argued in the paragraph following 
Eq. (P7D, a procedure involving the simple (RScl- and RSch-dependent) Borel transform 
leads in general to resummed predictions which can have significant dependence on the 
RScl and RSch used in the original TPS (including c 3 -dependence). Their approach (with 
r 3 = 130. and [2/2] dPA) would result in aJ^(Q 2 h ) = 0.2934jfto3?o for case I, and 0.2891±g;?gf| 
for case II; and af^(Mf) =0.1202±g;ggg for case I and 0.1140±g;^gf for case II. Comparing 
with results in Table |TI|, we see that these predictions are again very close to the predictions 
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of [2/2} 1 J 2 2 , the latter being based solely on the NNLO TPS ©. 

Recently, in the context of the Borel-Pade method of resummation (not used here), 
the knowledge of the location of renormalon poles was used in Ref. in two physical 
examples, to fix the denominator structure of the PA's of the Borel transform. 



V. BJPSR: USING PADE RESUMMED -FUNCTIONS 



Since nonperturbative physics appears to be of high relevance for the high-precision 
predictions in the case of the considered observable, one may go still one step further. Until 
now, we used for the /3-functions appearing in the integrated RGE ( JTJ| ) [cf. also (Q)] simply 
their TPS to the known order: 

TPS /3 (x) = -(3 x 2 (l + c\x + c 2 x 2 + c 3 x 3 ) , (59) 

where x = a s /Tc, and the bar over symbols denotes that they are different in different RSch's. 
However, in the nonperturbative region of large x, these TPS's may give wrong numerical 
results. To address this question, we may instead construct PA's based on these TPS's. PA's 
represent approximate analytic continuations (i.e., quasianalytic continuations) for the true 
f3(x) -functions from the perturbative (small-x) into the nonperturbative (large-x) region. 
A comprehensive source on mathematical properties of PA's is the book |45[. We have 
for fl5"9|) three PA candidates: [2/3]^, [3/2]^ and [4/1]^. Constructing these PA's on the 
basis of the TPS (j59f)) and then reexpanding in powers of x, gives us the higher order RSch- 
parameters Cj (j > 4) that were up until now simply set equal to zero. Only our approximant 
y/As^, and the ECH and TPS-PMS approximants for the NNLO TPS's (0), are sensitive 
to this change. Predictions a^ s (Qp h ) of Pade resummation approximants for S(Qp h ) in the 
previous Section, and the TPS evaluations themselves (NNLO, N 3 LO), are not affected by 
this change (they were calculated in MS RSch and at RScl Q 2 ] = Qph, nf — 3). 

For the approximant VAs 2 the relevant RSch's are those of a x (RSchl) and a 2 (RSch2), 
i.e., those with the RSch-parameters , c 3 , . . .) and (c| , c 3 , . . .), where the dots stand 
for Cfc and cj^ (k > 4) as determined by our choice of PA for the RSchl and RSch2 
/^-functions, respectively. Analogously, for the ECH and TPS-PMS approximants, the 
RSch's are (p 2 , C3, . . .) and (3p 2 /2, C3, . . .), where the dots stand for those RSch-parameters 
determined by our choice of the PA for the ECH and TPS-PMS /^-functions. So, each 
of the three choices of the PA defines, by the afore-mentioned mechanism of quasianalytic 
continuation into the nonperturbative sector, the unique schemes RSchl, ECH RSch, TPS- 
PMS RSch, and MS. 

For RSch2, we have to keep in mind one detail: In order to avoid presumably unnecessary 

(2) Ji) 



complications, the PMS conditions (|B.1|) -( [B~2"|) were written and used for the choice c 4 = c, 



1 



(5c4 = 0), so that the solutions fl39D-(f40D for Qi, Q 2 , dp and were independent of 
cz(=cP = c~P) and of all the other cjjf^ (k > 4; j = 1, 2). Therefore, once we choose a specific 
[M/N]p of the RSchl, the predicted c 4 must be reproduced also by the [M 1 /N'\p of the 
RSch2. This means that the order of the latter PA is by one unit higher than that of the 
former: M'+N' = M+N+1. Since the PA choices for the RSchl /^-function are [2/3], [3/2] 
and [4/1], those for the RSch2 /5-function are: [2/4], [3/3], [4/2], [5/1]. As to the numerics, 
the situation does not change much when different choices of [M'/N']p or even TPS for the 
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RSch2 are taken (with =c^\ and always the same fixed value of c 3 ). This is so because, 
in the strong-coupling regimes S > 0.155, ai is by a factor of 1.66 or more larger than a 2 . 
Concerning the choice of PA^ of MS RSch, this choice does not influence the predictions 
of c 3 at all, and influences only little the subsequent predictions for a!^ s (C}ph)- The latter 
is true mainly because of the hierarchy: a < a 2 < a\ (Q > Q 2 > Qi- Qi ~ 0.343C} p h) 
Q 2 ^0.672Q ph , Qo = Q ph «1.73 or 2.24 GeV). 

For the various PAp choices of RSchl, RSch2, ECH RSch and TPS-PMS RSch, we can 
just redo the entire calculation of the invariant Borel transforms Bs of ( |50"D and of their 
PA's, and find predictions for c 3 that give us the correct IRi pole y po \ e = 2. It turns out 
that the most stable c 3 -predictions in our approximant VAs 2 are those with [2/3)^1 for 
RSchl (01) and [2/4]^ for RSch2 {(32). The choice [2/3]^ and [5/1]^ gives virtually the 
same and almost as stable predictions for c 3 . For the ECH and TPS-PMS approximants, 
all three choices [2/3]^, [3/2]^, and [4/1]^ give comparably stable and mutually quite similar 
c 3 -predictions, but the choice [3/2] ^ seems to be slightly more stable than the other two. 
The results, for the mentioned optimal choices of PA^'s for the three approximants, are given 
in Table |TJ, in complete analogy with Table |. In some cases there are also other solutions 
for C3, not included in the Table, which differ significantly from those given in the Table. We 



PA^ 


c 3 for v^4: [2/%, [2/4]^ 


c 3 for ECH: [3/2]^ 


c 3 for TPS-PMS: [3/2]^ 


[2/1] 


21.7 


35.1 


35.1 


[3/1] 


15.7 


22.9 


21.5 


[4/1] 


15.8 


20.8 


18.7 


[5/1] 


16.9 


19.6 


17.3 


[1/2] 


12.8 


17.3 


17.3 


[2/2] 


14.9 


20.4 


19.4 


[3/2] 


15.8 


20.7±2.8i 


17.3±3.6i 


[4/2] 


15.7 


20.4±1.8i 


17.0±2.6i 


[1/3] 


15.0 


20.6 


19.5 


[2/3] 


15.1±1.2i 


19.3 


18.5 


[3/3] 


14.0±1.7i 


20.2±2.0i 


16.9±2.7i 


average 


ps 15.5 


~ 20. 


~ 19. 



TABLE III. As in Table 0, but the /3-functions in the approximants are taken as: [2/3] $ (RSchl), 
[2/4]/j (RSch2; cf^c^); [3/2]/j (ECH RSch, and TPS-PMS RSch). 



will adopt the approximate predictions as suggested by PA-g-'s of intermediate orders ([3/1], 
[4/1], [2/2], [3/2], [1/3], [2/3]): c 3 « 15.5 for y/A^; c 3 w 20 for the ECH; c 3 w 19 for the 
TPS-PMS. The actual values of c 3 must be exactly real. 

We recall that the results of the previous two Sections, including those of Table |, were 
for the simple choice of TPS/3 (|59"D for the corresponding RSch's ("truncated RSch's," with 
Cfc = for k > 4). Comparing those results with the results of Table |TJ, we see that the 
latter are somewhat higher and significantly more stable under the change of the choice 
of PA-g. This latter fact can be regarded as a numerical indication that it makes sense 



22 



to use certain PA resummations for the pertaining /^-functions of approximants when the 
considered observable (in this case BjPSR) contains nonperturbative effects. 

When the order of PA-g- is increased, the trend of the predictions is similar as in Table |: 
The predictions C3 tend to stabilize at intermediate orders of the PA-g's. The lowest order 
PA^'s ([1/2], and above all [2/1]) give unreliable predictions for c 3 , apparently because 
of a too simple structure of these PA's. The highest order PAg's ([5/1], [4/2], [3/3]) also 
sometimes give unreliable predictions, apparently because of their "overkill" capacity - these 
PA-g-'s depend on many terms in the power expansion of the approximant (up to ~S 7 ), while 
the original TPS (P^| ) on which the approximant is based is given only up to ~Oq (~a 3 ). 
Therefore, it seems plausible that the best and most stable predictions are given by PA^'s 
of intermediate orders ([3/1], [4/1], [2/2], [3/2], [1/3], [2/3]). 

With these choices for the values of C3 and for the pertaining /3-functions, we could now 
go on to calculating predictions of the three approximants for a^ s . Since the choice of PA/3 
for MS RSch will not matter much numerically, as we argued above, we could just choose 
blindly a PA^ or even the TPS for it. But at this point, we want to point out an additional 
argument for the made PAp choices of RSchl/RSch2, ECH RSch and TPS-PMS RSch. This 
argument will, in addition, lead us to a specific choice of PAp for MS RSch. 

In this context, we recall first that quasianalytic continuation, e.g. via PA's, of the TPS 
of a /3-function into the large-x (nonperturbative) region leads in general to a pole of such 
PAp(x) at some positive x. The authors of Ref. |46j pointed out that these poles "suggest 
the occurrence of dynamics in which both a strong and an asymptotically-free phase share a 
common infrared attractor." Now, if there is such a common point x po \ e = a^ ole /tt where the 
two phases meet, it is reasonable to expect that its numerical value does not vary wildly when 
we change RSch - provided that the RSch's in question are themselves physically motivated 
(physically reasonable) in the nonperturbative regime.^ Such physically motivated RSch's 
should include those connected in some significant way with the calculation of the considered 
observable and of the predicted coupling parameters. In the case of our approximant y/A<p, 
these are RSchl and RSch2, and in addition MS when we want to extract a^ s (Q p h) from 
the approximant. In Fig. |2| we present the TPS's of RSchl, RSch2 and MS /3-functions, as 
well as the previously chosen [2/3]^ of RSchl and [2/4]^ of RSch2 (cf. Table |T|; C3 = 15.5), 
and we include also [2/3]^ of MS RSch. The Figure shows that all these PA /3-functions 
have about the same x po i e (x po i e = 0.334, 0.325, 0.311, respectively). The mutual proximity 
of Xp ie's of RSchl and RSch2 PA^'s is now yet another indication that these PA^'s, chosen 
previously on the basis of the stability of C3-predictions, are the reasonable ones. Further, 
[2/3]/? appears to be the reasonable choice for MS RSch. The choices [3/2] ^ and [4/1]^ for 
MS RSch give x po i e = 0.119, 0.213, respectively, which is further away from the x po i e of RSchl 
and RSch2. We could choose, in principle, for RSchl and RSch2 other PA^'s. We recall that 
for RSchl we can have: [2/%, [3/2]^, [4/%; for RSch2: [2/4]^, [3/3]^, [4/2]^, [5/1]^. 
However, when taking [3/2] or [4/1]^, we always end up either with a situation when 
the two positive x po ie values of (31 and (32 are far apart, or both are unphysically small, 
or one positive x po i c doesn't exist, or there are virtually no predictions for C3 (not even 
unstable ones), or x pole values are very unstable under the change of c 3 in the interesting 



In the perturbative regime, all RSch's are formally equivalent. 
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region c 3 ~ 12-16. Concerning the latter point - when taking [3/2]^, and for RSch2 [3/3]^ 
or [4/2]^, the location of x po i e of the latter PA^'s changes drastically when c 3 is varied 
around the interesting values of 12-16, thus signalling instability of these PA^'s. The choice 
[2/3] 0i and [5/11/32, which gave very similar and almost as stable results for c 3 as the most 
preferred choice [2/3] ,31 and [2/4] 02 , gives the corresponding poles again close to each other: 
x po ie = 0.334, 0.291, respectively. So, the PA^ choices [2/3]^ and [2/4] ,32 (or [5/1]^) f° r our 
approximant give us the most stable C3-predictions and are the only ones giving mutually 
similar (and reasonable) values of x po i e of RSchl and RSch2. 

It is also encouraging that the choices [3/2] ,3 for the ECH and TPS-PMS RSch's give 
us Xp i e values comparable to the ones previously mentioned: a; po i c = 0.263 for ECH with 
c 3 = 20; x po i e = 0.327 for TPS-PMS with c 3 = 19. Even other choices of PA^ for the ECH and 
TPS-PMS RSch's ([2/3]/3, [4/1]^), which also gave rather stable and similar C3-predictions, 
give us Xp i e ~ 0.27-0.41. Hence, also in this case we see correlation between the stability 
of the c 3 -predictions on the one hand and x po i e ~ 0.3-0.4 on the other hand. 

The authors of Refs. f4"T|j4"8|j estimated the 5-loop coefficient c^ s of the MS /3-f unction, 



by applying their method of Asymptotic Pade Approximation (APAP, j|7|) and its im- 
provement using estimators over negative numbers of flavors (WAPAP, ||48|| ). Their pre- 
dicted values by two variants of the latter method, when including the four-loop quartic 
Casimir contributions, are c^ s = 123.7, 115.3 (cf. Tables III and IV in Ref. |4"8|] , respectively; 



n/ = 3). On the other hand, the simple PA's [2/3], [3/2], [4/1] for MS /^-function predict 
= 6 2.2, 149.8, 98.5, and x po i e = 0.311, 0.119, 0.213, respectively. If we assume that the 



actual value of c 4 is close to the one predicted by [48], and if we were led just by the 



requirement that the PA should reproduce well this value, then [4/1] would be the preferred 
choice. However, the authors of j|8| indicated that their predicted value of c 4 may be changed 
significantly if new Casimir terms, appearing for the first time at the 5-loop order, are large. 
Our choice [2/3] for MS /3-function was motivated by the value of x po ie = 0.311 lying close 
to Xpoie of the /3-functions appearing in the discussed approximants for the BjPSR. Further, 
the precise choice of the PA for MS /3-function practically does not influence the numeri- 
cal results of our analysis, because ao = a(lnQp h ; df , . . .) is significantly smaller than the 

coupling parameters aj = a(\n Q 2 -\ c 2 J , c 3 , c 4 , c 5 , . . .) (j = 1, 2) appearing in our approximant, 
and the parameters oech and apMs appearing in the ECH and the TPS-PMS approximants. 
To summarize: 

• the best choice in calculating a^ s from our approximant \CA<p is: c 3 ~15.5; the PAp 
choice [2/3]^ for RSchl, [2/4]^ choice for RSch2 {cf =c£ ) ); and [2/3]^ for MS RSch; 

• the best choice in calculating from the ECH and TPS-PMS approximants is: 
c 3 w 20 and 19, respectively; the PAp choice [3/2]^ for ECH RSch and TPS-PMS 
RSch; and [2/3],? for MS RSch; 

• our, the ECH and the TPS-PMS approximants are completely independent of the 
original choice of the RScl and RSch, because the c 3 parameter is determined by using 
the RScl- and RSch-invariant Borel transform B(z) of Sec. III. 



In practice, this means that for our approximant -^/As 2 the two coupling parameters 
,2. 
j- 



dj = a(lnQ|; c£ \ c 3 , c 4 , Cg , . . .) (j = 1,2) are now related with the coupling parameter 
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MS JAS „MS 



a = a(\nQ 2 ;cf s ,c 
equation gl]) [cf. also Q -Q]: 



via the following (PA-)version of the subtracted Stevenson 



A)ln(^) =- + Cl ln 



+ 



{PA^x) + (3 x 2 (l+ Cl x)} 



o 



;i + cix)PA /3i (x) 



1 



Ci In 

a 



Cl a 



l + cia J 



°o^ {[2/3] MS /3(x) + f3 x 2 (l + Cl x)} 
o x 2 (l+cix)[2/3] M s/3(a;) 



(60) 



where PA^- stands for the mentioned [2/3]^ of RSchl (when j = 1) and [2/4]^ of RSch2 (when 
j = 2), with C3 = 15.5. We recall that the scales Q 2 and the parameters c 2 (j = 1,2) of 



the approximant, which are RScl- and RSch-invariant and calculated in Sections II and III 
[cf. (|3T?|)-(|4"0"|)], are independent of the parameter c 3 and of any higher order /3-parameter 
(k > 4; j = 1, 2) appearing in aj =a(\nQ 2 - ci\ c 3 , c 4 , cf, . . .). For the ECH and TPS-PMS 
the calculation is performed in an analogous way. 

The results of these calculations, i.e., the predicted values of a^ s (Qph) an d a^ s (M§), 
are given in Table |Iy| for the approximants VAs 2 ) ECH and TPS-PMS. The predictions are 



approximant (with PA^'s) 


a s (5 GeV 2 ): (I) a s (3 GeV 2 ): (II) 


a s (M 2 ): (I) a s (M 2 ): (II) 


\/A% 2] (c 3 = 15.5) 


fl ooqo+0.0182 n 98n c+0.0297 


n 1 1 QV+0.0028 n 1 9 7 +0.0041 
-0.0054 U - ii2 '-0.0189 


ECH (c 3 = 20.) 


n 9S^R+ U - U195 n 90O9+U.U325 
U.Z<3OO_ 0321 U.ZOZZ_ 0993 


u.iiyu_ 00056 u.ii3U_ 00192 


TPS-PMS (c 3 = 19.) 


0.2867± ; 0328 0.2831 ± 1001 


n 1 09+0.0030 n 1 1 01 +? 

u - iiyz -0.0057 u ' lioi -0.0192 



TABLE IV. Predictions for for our, ECH, and TPS-PMS approximants, when the PA-resummed 



/^-functions in the approximants are taken as in Table III. Predictions for the case I (|53|) and II (|5^) are 
given in parallel. 

now a little, but still significantly, lower than those of the corresponding approximants in 
Table |T| where all the /^-functions were taken in the TPS form ( |59D and with c 3 = 12.5, 17, 16, 
respectively. The evolution from a^^Qph) to a^ s (M|) was performed as in the previous 
Section, i.e., with the four-loop RGE (i.e., TPS /3-function of MS) and the corresponding 
three-loop flavor threshold matching conditions. If we replace the TPS /3-function of MS by 
its PA [2/3]^ in the RGE for the evolution af^iQ 2 ^) i-> d^(M%), the results for ctp(Mf) 



decrease insignificantly (by less than 0.04%) and the numbers in Table [IV| do not change 



In Fig. |3|(a) we present predictions S(Q 2 h ) as functions of a^ s (Qph) ( n / = 3, e.g. Q 2 h = 
3 or 5 GeV 2 ), and in Fig. |(b) the predictions for S(5GeV 2 ) as functions of crp(Mf), 
for the three approximants with the afore-mentioned PA choices for the /9-functions. For 
comparison, we include in the Figures also predictions of these three approximants when 
all the /9-functions have the TPS form (|59"D and the correspondingly smaller C3's (the latter 
curves are contained also in Figs. [I]). Predictions of the PA resummation approximants (for 
S) are not included, since these methods are insensitive to the mentioned PA-quasianalytic 
continuation of the /^-functions and the results remain for them the same as in Figs. [I] 
and Table |H]. We presented in Figs. § the curves for the case of approximants with the 
mentioned PA /3-functions only so far as the method works. More specifically, when the 
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integration interval in the first integral of (pOD starts including values x larger than those 
at which the absolute value of the PA^ exceeds the value 2, we stop the calculation of the 
approximant since the latter would otherwise probe values too near the pole of PA^ (i.e., 
too near the common point of the asymptotically-free and the strong phase) and would thus 
be unreliable. 

The considered BjPSR observable S(Qp h ) has a higher-twist (h.t.) contribution, esti- 
mated from QCD sum rule ||49[|p| 

s<-wy « (a09±0 n ° 45)GeV2 . («) 

Vph 

which should be added to the perturbation series for S. If adding this term in the numerical 
analysis, the predicted central values of a^ IS (M|) given in Table [Ij decrease significantly. 
For example, the NNLO TPS central value predictions (M§) = 0.1252 (case I) and 
0.1183 (case II) then decrease to 0.1200-0.1236 (case I) and 0.1091-0.1157 (case II), where 
the lower and upper values in each case correspond to the largest and the smallest value 
choice in (|H|). This indicates numerically that our approximant (c^ = 15.5, Table fV\j , which 
gives the central values af s (M 2 z ) = 0.1187 (case I) and 0.1127 (case II), already contains at 
least part of the nonperturbative effects from the leading higher-twist operator (~1/Qp h ). 
The same is true for the ECH (03 = 20.) and TPS-PMS(c3 = 19.). In order to understand 
this numerical indication, we recall that the information on the location of the leading IR 
renormalon (IRi) pole of the considered observable has already been incorporated in these 
approximants, via the afore-mentioned fixing of the value of C3-parameter. And the so called 
ambiguity of the leading IR renormalon is of the same form ~1/Qp h as the higher-twist term 
(|Bl"l), and even the estimated coefficients are of the same order of magnitude |5l] (cf. also 
Ref. p6[ on this point). Our approximant, the ECH and the TPS-PMS, via the discussed 



C3-fixing, implicitly provide approximant-specific prescriptions of how to integrate in the 
Borel integral over the IRi pole, thus eliminating the (leading) renormalon ambiguity. 

VI. DISCUSSION OF THE NUMERICAL RESULTS 

The main reason to apply our approach (and PA approaches) to the BjPSR was to 
investigate efficiencies of various methods and the influence of the nonperturbative sector. 
Another reason was that the BjPSR is a Euclidean observable (q 2 h = —Qp^ < 0), and for such 
observables various resummation methods are believed to work well since no real particle 



thresholds are involved in the observable 52.53 



The main prediction of our approximant V 'As* can be read off from Table fVj, for two 

2 

P h 



cases (f)3|) and (ET) of the BjPSR-integral values at Q^ h = 5 and 3 GeV 2 , respectively. 



extracted from experiments 

aT\Ml) = 0.11871™ (/); 0.1127^ (H). (62) 



13 Deficiencies of the QCD sum rule calculations were pointed out in 
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The ECH and the TPS-PMS give results similar to these, when C3-parameter in them is 
adjusted in the afore-mentioned way - see Table IV. The diagonal PA (dPA) methods give 
higher predictions, and the nondiagonal PA methods even higher - see Table [TI] and Figs. [I]. 

The result (^) for case II, which is based on the measurements before 1997 and a NLO 
pQCD extrapolation for low x-q] |37| (|52|), shows quite large uncertainties, a consequence 
of the large uncertainties (p6|) . (]52])1. The result (|^) for case I, based on the most 
recent measurements and a similar NLO pQCD extrapolation for small XBj, by the SLAC 



E155 Collaboration |19j (|53|), already shows significantly reduced uncertainties. This is so 



to a large degree because of additional new measurements in the Iow-xbj regime. And most 
importantly, the central values of case I in (|62| ) are now significantly higher than those of (the 
older) case II. We recall that the central values in ( |62|) correspond to the central values of the 
BjPSR-integral (|5tJ) and ([52]). We did not attempt to estimate the theoretical uncertainties 
originating from the resummation method itself. However, the combined results of Table [IV| 
a^ s (M|) = 0.119to;oog f° r ( new ) case 1 could be regarded as containing nonconservatively 

estimated theoretical uncertainties. 

The present world average is o:f s (Mf) = 0.1173 ± 0.0020 by Ref. 0], and 0.1184 ± 
0.0031 by Ref. ||55|| . Predictions of the simple (NNLO) TPS evaluation in (new) case I give 

0. 1252^0 0078 (see Table [TJ), the central value and most of the interval lying significantly above 
the world average. On the other hand, the simple (NNLO) TPS evaluation in (older) case II 
predicts 0.1183io;o232 ( see Table [IT]), the central value agreeing well with the world average, 
but the uncertainty interval being much broader. However, the situation changes drastically 
when employing more sophisticated resummation methods. The values for BjPSR-predicted 
a^ s (M|) go down the more significantly, the more sophisticated resummation we perform - 
cf. Table [H]for the PA-methods, and for TPS-PMS, ECH and y/As^ when the /3-functions 
have truncated form, and Table [TV] for the last three methods when the /3-functions are 
resummed. The predictions of approximants in the latter Table have a^ s (M|) « O.lWj^oog 
(case I, new) and 0.113^019 (case II, old). The predictions of (new) case I now agree well 
with the world average 0.1184 ± 0.0031 of Ref. ||55|| , while those of (older) case II lie almost 
entirely below the world average intervals. 

Thus, the use of resummation methods which account for nonperturbative contributions 
by the mechanism of quasianalytic continuation and by incorporation of the information on 
the leading IR renormalon pole, predict the values of a^ s (M|) which agree well with the 
present world average if the most recent BjPSR data |19| are used. This suggests, among 
other things, that for reliable predictions of a^ s from reasonably well measured low-energy 
QCD observables, we have to know the NNLO terms (~a 3 ), employ nontrivial resummation 
methods, and possibly incorporate some nonperturbative (renormalon) information in the 
resummation. 

Some of the recently performed analyses beyond the NLO, by other authors, gave predic- 
tions: crp(Mf) = 0.118 ±0.006 fl5§] from the CCFR data for x Bi F 3 structure function from 
WVDIS (NNLO); 0.1172±0.0024 [§7|] from £N DIS (NNLO); 0.115±0.008 |58| and 0.114±g;g^ 
f59| from Gross-Llewellyn-Smith sum rule (NNLO); 0.1181 ±0.0031 from hadronic r-decay 
(NNLO, combined results, [0); 0.115 ±0.004 p0|j54]] from lattice computations. 

We note that the BjPSR predictions deviate from the world average in case I' fl54"l), 

1. e., when we include in the experimental data of case I the nuclear effects originating 



from spin-one isosinglet 6-quark clusters in deuteron and helium according to Ref. [41 
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on top of the nuclear wavefunction effects and NLO pQCD small-XBj extrapolation effects: 
ap(M|) w 0.103l8;^7 (NNLO TPS); 0.10ll£[j| (dPA, ECH, TPS-PMS, our approximant). 
The combination of (older) case II results and the mentioned 6-quark cluster nuclear effects 
(case II') increases the value of the BjPSR integral so much that the predicted values of 
a^ s (M|) are unacceptably low: the central values would be 0.094-0.095 for all approxi- 



mants; the maximal allowed values would be about 0.113 by the methods of Table [TV] and 
0.114 by the dPA. 



The authors of Ref. |37] obtained, among other things, the BjPSR-predicted values 

afS(M§) = 0.1181^26, apparently using the simple NNLO TPS sum (g§) directly in their 

analysis. They used the BjPSR-integral values (|52|) , i.e. here case II, which were extracted 

by them from low-Qp h SLAC experiments carried out before 1997. They used the value of 

1 5^1 = 1.257 known at the time, in contrast to the value of (|38|) . Their RGE evolution from 
Q2 h = 3GeV 2 toM | 

was apparently carried out at the three-loop level, since the fourth-loop 



/^-coefficient c 3 (rif) [H2] and the corresponding three-loop flavor-threshold matching |53[ 
were not known at the time. These two effects largely neutralize each other and their result 
is then close to the NNLO TPS result for case II (Table 0): af s (Ml) = 0.118 +0010 



0.023- 

The authors of Ref. |6| obtained the BjPSR-predicted values af S (M|) = 0.1161^05 ± 
0.003. They used a dPA method of resummation [2/2] s mentioned towards the end of Section 
IV. However, they took the BjPSR-integral values ([51]) where the naive Regge small-a^j 
extrapolation was used, and apparently the value \g&\ = 1.257 known at the time. Further, 
they included the effects of the higher-twist term ( |61| ) on top of their dPA resummation. 
The additional uncertainty ±0.003 can be called the method uncertainty. It was estimated 
by them by additionally using the results of the nondiagonal PA resummations [1/2] 5 and 
[2/l]s, the RScl-dependence of their dPA results, and the uncertainty of the higher-twist 
term. 

When we reexpand the approximants in powers of the original ao (at RScl Qq = Qp^, in MS 
RSch, rif = 3), we obtain predictions for coefficient r 3 at <2q of expansion fl3~6|) - cf. Eq. ([47]) and 
the discussion following it. Our approximant, with c 3 = 15.5, predicts r 3 = 125.8— c| IS /2+C3 « 
130.8. The ECH approximant, with c 3 = 20., predicts r 3 = 129.9 + (-cp+Cs)/2 w 129 .4 
The two predictions are close to each other, suggesting r 3 = 130. ± 1. This agrees well with 
the prediction of Ref. |52j] r 3 ^ 129.9 (~130.) which was obtained from the ECH under the 
assumption (— dj +C3) ~0 (note that c 3 vIS ~21.0 [42| was not even known at the time Ref. |52| 
was written). 

The predictions for r 3 , as well as the values of Q\, Q%, c 2 , c| (p9|)-(^0|) and of c 3 (Tables 
|, PT| ), are for rif = 3 and are, of course, independent of the specific values for the BjPSR 
integral (|53"D, (|52|) [<=>■ (^5|)-(0)] that we subsequently used to obtain values for a^ s (M|). 



VII. SUMMARY AND OUTLOOK 



We presented an extension of our previous method of resummation JI5|^TT|| for truncated 



perturbation series (TPS) of massless QCD observables given at the next-to-next-to-leading 
order (NNLO). While the previous method, partly related to the method of the diagonal 
Pade approximants (dPA's), completely eliminated the unphysical dependence of the sum 
on the renormalization scale (RScl), the extension presented here eliminates in addition 
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the unphysical dependence on the renormalization scheme (RSch). The dependence on 
the leading RSch-parameter = 02 / flo is eliminated by a variant of the method of 
the principle of minimal sensitivity (PMS). The dependence on the next-to-leading RSch- 
parameter cf^ = (3^ I (3q is eliminated by fixing the C3-value in the approximant so that the 
correct value of the location of the leading infrared renormalon (IRi) pole is obtained (by 
PA's of an RScl- and RSch-invariant Borel transform). Hence, in the approximant we use 
/^-functions which go beyond the highest calculated order in the observable (NNLO) - in 
order to incorporate an important piece of nonperturbative information (IRi pole location) 
which is not contained in the available NNLO TPS anyway. The results are apparently 
further improved when we resum those /3-functions which are relevant for the calculation 
of the approximant (RSchl and RSch2 /3-functions, for a\ and a 2 ) and of c^ s {Qph) (MS 
RSch), by judiciously choosing certain PA-forms for those /3-functions. 

We applied this method to the Bjorken polarized sum rule (BjPSR) at low values of the 
momentum transfer of the virtual photon Q^ h = 5 or 3 GeV 2 . The C3-fixing by the IRi pole 
location is well motivated in this case, because the contributions of the leading ultraviolet 
renormalon (UVi) appear to be sufficiently suppressed in comparison to those of the IRi. 
We compared predictions of our resummation with the values for the BjPSR integral (|53|) 
and (|52"D extracted from experiments, and obtained a^ s (M|) = 0.1187lo!oo54 (new case I) 
and 0.1127lo;o?gg (older case II), respectively. Here, the central values 0.1187 and 0.1127 
correspond to the central values in (^) and (0), respectively. For more discussion on the 
issue of the experimental values (|53|) and (|2|) (cases I, II) we refer to Sections IV and VI. 
It is gratifying that the newest available experimental values ([53]) lead to predictions for 
which agree well with the present world average. The results of Grunberg's method of 
the effective charge (ECH) and of Stevenson's TPS-PMS method give very similar results 
(cf. Table |P7D if the C3-parameter in these methods is fixed by the same afore-mentioned 
requirement as in our approximant and PA-forms of the pertaining /5-functions are chosen 



analogously. The combined result of Table IV, in case I, i.e. with the newest data of Ref. |19 
is 

of (Ml) = . (63) 

The dPA methods of resummation of 5* predict higher values (central values about 0.120 in 
case I; 0.114 in case II), the non-diagonal PA's even higher (central values about 0.122 in 
case I; 0.115 in case II), and the NNLO TPS itself the highest values (central value about 
0.125 in case I; 0.118 in case II). 

We expect that our approximant \CA<p, as well as the ECH and TPS-PMS, produced 
reliable resummation results for the considered observable, because - via their dependence 
on C3 - we can incorporate into them in the afore-mentioned way important nonperturbative 
information about the IRi pole, and simultaneously achieve full RSch-independence. The 
C3-dependence in \/As%, in the ECH and in the TPS-PMS, is very closely related with the 
sensitivity of these approximants to the details of the corresponding RGE evolution. These 
details (c3-terms) in the RGE evolution are numerically more important in the lower energy 
regions, i.e., when the relevant energies for the observable are low. Thus, significant C3- 
dependence of these approximants signals the relevance of nonperturbative regimes for the 
observable [cf. Eqs. (|32})-(|3~3|)1. It then appears natural that the c 3 -parameter in these 
approximants, i.e. the only parameter left free, is made to parametrize the location of the 



29 



(nonperturbative) IRi pole. The (d)PA's, in contrast, possess besides the C3-dependence also 
dependence on the leading RSch-parameter C2, and even on the RScl. Thus, the parameter 
C3 in them is not in a special position, and there is more ambiguity as to how to incorporate 
into the PA's the information about the IR X pole. 

It appears that the leading higher-twist term contribution to the BjPSR (~l/<5p h ), 
or a part of it, is implicitly contained in VAs 2 ; as well as in the ECH and the TPS- 
PMS, via the afore-mentioned C3-fixing. In this context, we point out that the so called 
renormalon ambiguity arising from the IRi of the BjPSR has the form ~1/Qph, i.e., the 
form of the leading higher-twist term. Even the coefficients of this term, as estimated by 
the renormalon ambiguity arguments, are of the same order of magnitude as those predicted 
(estimated) from QCD sum rule. One can say that the described approaches implicitly 
give approximate-specific prescriptions for the elimination of the (leading IR) renormalon 
ambiguity. 

Looking beyond the numerical analysis of the BjPSR, we wish to stress that in cases of 
other QCD observables that are (or eventually will be) known to the NNLO, the analogous 
numerical analyses may give different hierarchies of numerical results. Actual resummation 
analyses should be performed also for such observables, in order to shed more light on the 
questions about the relative importance of various kinds of contributions. 

The (d)PA methods, when applied directly to the (NNLO) TPS's, are trying to include 
some nonperturbative contributions through quasianalytic continuation of the TPS from 
the perturbative (small-a) to the nonperturbative (large-a) region. In the course of this 
continuation, the pole structure of the Borel transform of the sum may be missed, but some 
other nonperturbative (but less singular) features of the sum itself may be reproduced well. 
But our approximant s/As* would presumably do at least as good a job as the dPA's in 
reproducing these latter nonperturbative features. This is so because \/A~s^ (|L|) reduces 
to the dPA [2/2] % 2 in the large-/?o (one-loop RGE evolution) approximation when thus 
the full RScl- and RSch-invariance requirements are abandoned - cf. discussion following 
Eqs. ©-([□]). The ECH and the TPS-PMS methods do not possess this strong "[2/2] 1 / 2 
PA-type" mechanism of quasianalytic continuation, since these two methods fix the RScl 
and the RSch in the TPS itself without going beyond the (NNLO) polynomial TPS form 
in a. The ECH, and somewhat less explicitly the TPS-PMS, possess a weaker type of 
quasianalytic continuation, because the one-loop RGE-evolved a = a s /ir (from ao) is a [1/1] 
PA of a . 

Stated differently, our (NNLO) approximants, from a theoretical viewpoint, combine the 
favorable feature of the (d)PA's (strong quasianalytic continuation into the large-a regime) 
with the favorable feature of the TPS-form NNLO approximants ECH and TPS-PMS (full 
RScl- and c 2 -independence). The residual RSch-dependence (c 3 -dependence) in the latter 
approximants and in our approximant allows us to incorporate into them, often in a well- 
motivated manner, nonperturbative information on the location of the leading IR renormalon 
pole, and to achieve in this way simultaneously the full RSch-independence as well. 
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nuclear effects in the BjPSR. 




Appendix A. EXPANSION OF THE GENERAL COUPLING a IN POWERS OF a 

We outline here the derivation of the expansion of QCD coupling a = a(\nQ 2 ; c 2 , C3, . . .) 
(a = a s /n) in power series of a = a(lnQl; c 2 \ cf\ . . .). The starting point is the Stevenson 
equation (cf. Ref. ||, first entry, Appendix A) which is obtained by integrating RGE (^) 

1 . / cio \ f a , \ 1 1 
- + ci In + / dx -r— - - : 

a W+cidJ Jo [x z {l + cix) x z {l + cix+C2X z + csX' i + - ■ ■ 

(A.l) 

It can be shown that A here is a universal scale (~0.1 GeV) independent of the scale Q and 
of the scheme parameters Cj (j >2). Writing the analogous equation for a , and subtracting 
the two, we obtain 

/3 ln ( 7L ) = - + c\ In ( - Cl<X \ + / dx— ( c 2 + c 3^ + ) 

\QoJ a Vl+cio/ J {l+c 1 x){l+c 1 x+C2X z +c 3 x' i -\ ) 

^ - * m (^) - r* ( ^ + "°, m .. ■ • (A.2) 



a Vl+cia y Jo (i + ClX )(i +ClX + c (°) x 2 + c (o) x 3 + 

This equation determines a as function of clq. The solution a in form of a power series of ao 
is the Taylor series for function a of multiple arguments In Q 2 and c/s (j > 2). To obtain this 
power series, one way would be to find first the derivatives da/dcj [the derivative da/dlnQ 2 
is already given by RGE (^)j. For this, we take the partial derivative of both sides of the 
above equation with respect to Cj (j > 2) and after some algebra we obtain 

— = a(l + c 1 a + c 2 a + c 3 a + ■ ■ •) J q (1 + ~ + c ^ + ^3 + . . . )2 • (A-3) 

Expanding the integrand in powers of x and integrating out each term, we obtain the partial 
derivatives as power series 

*-rf(l + S* + ...). (A.4, 

■)■ < A - 5 > 

Repeated application of these equations, as well as of RGE (|) itself, leads us to the following 
Taylor expansion of a in powers of a = a(\nQl; Cg , cf\ . . .): 

a = a + Oq(— x) + ag(x 2 — C\X + 5c 2 ) 

+ag(-a; 3 + ^cix 2 - c^x - 3x6c 2 + 7^3) 



x 4 -^c 1 x 3 + (^c 2 + 34 0) + 6 ( 5c 2 )x 2 



+(-4 0) - 3 Cl Sc 2 - 25c 3 )x + {\cf8c 2 + |(<5c 2 ) 2 - i Cl 5c 3 + he,)} + C?(ag) , (A.7) 

3 3 O 3 J 
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where we denoted 



a(lnQ 2 ;c 2 ,c 3 , 



a = a {lnQ 2 ]c { 2 \cf\. . .) 



Q 2 

x = Pom—*, 
Qo 



5c k = c k - cf> 



(A.8) 
(A.9) 



Appendix B. EXPLICIT PMS CONDITIONS 

Here we will write explicitly the PMS-like conditions (|31~1) in its lowest order 



To do this, we calculate explicitly the derivatives (^Tj) and then expand them in powers of 
d = a(lnQl; c 2 = c 2 ^; C3; . . .) to their lowest nontrivial order .0 We assume relation (|34"D, 
i.e., 5c3 = 0, and in addition 5c^{= C4 — c&) = 0. Further, we use relations (|2~6|)-(p7|) and 



notations (|28|) -(p0|). The results, obtained with help of Mathematica, are the following: 



OA 



[2/2] 



dd s) 



= -d b J 27(<5c 2 ) a - 157 Cl (5c 2 ) V - 85c 2 yt -27ci + 12c { 2 s) + 34yi-8z 2 (c 2 S) ) 



.2 1' Oh 



+4Sc iy d _ 13y 2 _-3z 2 (c { 2 s) ) \\6yi 5 Cl ck 2 + 16^ -8^(c 2 s V) \ +O(d b ) = 



(B.l) 



54 



[2/2] 



«9(5c 2 ^ 



'2 



27(5c 2 ) 4 - 315 Cl (5c 2 ) V + B^c^)^ 7c 2 -2cf +3^(4' 



00' 



-80 Cl 5c 2 y_ [-24 s) ^-24 s) z 2 (c 2 s) ) + 12^(4 s) ) 2 /! + 34(c 2 s) ) + 7c? (jd+zg(c 
+ 12(5c 2 ) 2 [-24V-24 s) ^(4 s) ) + 15z 2 (c 2 s) )y 2 +3^(4 s) ) + c? (82y 2 +7^)) 



X<{12^ 5c 1( 5c 2 + 16^-8z 2 (c 2 s) )i/ 



0{a\ 



0) 







(B.2) 



The actual PMS-type equations are now obtained by requiring that the coefficients at ~a[j 
in ( |B.l| )-( [Tj~2"| ) be zero. When we have several possible solutions of the coupled system ( |26| ) 
and (|B.1|)-(|B~2|) for the three unknowns y_, c 2 and 5c 2 , we have to choose, in the PMS- 



spirit, among the resulting approximants that one which has the smallest curvature. The 
curvature can be calculated by first obtaining the eigenvalues CAi and CA 2 of the curvature 
matrix Ca'- 



a 2 A s 



d 2 A e 



a 2 A? 



a 2 As 



a4W a(4 2) )2 



(B.3) 



'CAi 
,CA 2j 



1 d 2 As 



+ 



d 2 Ao 



4«9(c 2 s) ) 2 d(Sc 2 ) 



± 



( d 2 A- s 
\d(5c 2 )dc 



(«) 



+ 



1 d 2 As 
4«9(c 2 s) ) 2 



d 2 A § 
d(5c 2 y 



1/2 



(B.4) 



14 In fact, a with any RScl and any RSch-parameters would do the job and give the same coefficient 
at the leading nontrivial order a 5 . 
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In the last expression, we traded c 2 and c 2 for c 2 = (c 2 + c 2 )/2 and 5c 2 = (c 2 — c 2 ). 

[2 /2l 

The curvature C4 of the solution A S L can be defined in at least two obvious ways which 
are virtually equivalent 

C A = \CA 1 \ + \CA 2 \ , or: C A = yfiCAtf + {CA 2 f . (B.5) 



Appendix C. ECH AND TPS-PMS METHODS FOR NNLO TPS 

The effective charge method (ECH) || of resummation of the NNLO TPS S\ 2 ] (0) can 
be expressed by employment of the subtracted version ( |A.2| ) of Stevenson equation 



1 . / cia \ f a ° , (4 0) +4 0) x + ■ ■ 

-ri + — + a In - + / dx- ' 



ao Vl+cia / Jo (l + cix)(l + cix+4 0) x 2 + 4 0) x 3 + • ■ ■) 

Cl in ( 7 iqech ) + r cH dx 7 — (pi+^+---) — . (ai) 



ciech Vl + cia E cH/ -/o (l + cia;)(l+cia;+p2a; 2 + C32; 3 + ■ 

The ECH resummation value is S' ECH = oech- In superscript "(0)" denotes the 

original RSch of S[ 2 ] (for example MS RSch with rif — 3 in the considered BjPSR case), and 
C3 denotes the NNLO ECH value of C3 (in principle unknown at NNLO). Further, c 2 = p 2 , 
the latter RScl- and RSch-invariant is defined in ( p^f) . The coupling ao^a^/Tr is defined 
a = a(ln(5oi c 2°\ Cs ^) • • •) as m (!); Qo being the original RScl in the TPS (chosen equal 
3 GeV 2 in the considered BjPSR case); t\ = — ft ^(Qeck/Qo) ^ s * n e NLO TPS coefficient 
as staying in (Q) at the original RScl Qq. In the above relation ( |C.1| ), we often ignore the 
terms oc and Ck (k > 4) since they are not known, i.e. we often choose the TPS form 
for the /3(x)-functions. For a given value of do, solving the above relation numerically for 
Qech gives us the resummed prediction for observable S. It is dependent on C3 which, at 
this stage, is not known. More explicitly: 

S ECn (c 3 ) = a E cu(c 3 ) = a(\nQl CH ;p 2 ,c 3 , . . .) , with : Q| CH = Q 2 exp(-ri//3 ) • (C.2) 



For the TPS-PMS method § applied to the NNLO TPS S [2] , relation still remains 

valid, but with the replacements 

a EC H(c 3 ) ^ a P Ms(c 3 ) , c ECH = p 2 h-> c^ MS = ^p 2 . (C.3) 
The resummed expression in the (NNLO) TPS PMS case is the following TPS: 

S PMS (c 3 ) = apMs - ^P2o|ms ' with : a PMs(c 3 ) = a(lnQ| CH ; (3/2)p 2 , c 3 , • • •) , (C.4) 

which again depends on c 3 . Expression ( |C.4| ) is obtained by imposing PMS conditions on 
the TPS S [2] (\nQ 2 ; c 2 , c 3 , . . .) = S PMS : dS [2] /d\nQ 2 ~ a 5 ~ dS [2] /dc 2 . It is straightforward 
to verify that, if p 2 > (as in the considered BjPSR case), S FMS is bounded from above due 
to its specific TPS form: 5* PMS < (2/3) 3 ^ 2 p 2 1 ^ 2 , which in the considered BjPSR case (|36|) is 
0.2326 (because p 2 = 5.476). 
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T 




a s (M|) 

FIG. 1. Predictions of various approximants: (a) for S(Qp h ) as functions of Q^ s ((5p h ) when 
nj = 3; (b) for S(Q^ h = 5GeV 2 ) as functions of a^ s (M|). The values of the C3-parameter in our 
approximant (03 = 12.5), ECH (03 = 17.) and TPS-PMS (03 = 16.) have been adjusted to ensure 
the correct location of the leading IR renormalon pole. The experimental bounds S m i n , S maiX and 
'S'mid are indicated as dashed horizontal lines for case I (||) (Qp h = 5GeV 2 ) and dotted horizontal 
lines for case II © (Q 2 h = 3GeV 2 ). 
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FIG. 3. Predictions of our approximant (with c 3 = 15.5), ECH (with c 3 = 20.), and TPS-PMS (with 
c 3 = 19.): (a) for S(Q 2 ph ) as functions of af^iQ 2 ^) when n f = 3; (b) for S{Q 2 ph = 5GeV 2 ) as functions of 
a^ IS (Af|). The PA choices of the RGE /3-functions were made as explained in the text. For comparison, we 
include also the corresponding predictions from Figs, [l] when the TPS's (|5^) are used for the /3-functions. 
The values of the C3-parameter have been adjusted in all cases to ensure the correct location of the leading 
IR renormalon pole. The experimental bounds are denoted as in Figs. m. 
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